Electroproduction of r\ Mesons in the 5n(1535) Resonance Region at High Momentum 
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The differential cross-section for the process p(e, e'p)n has been measured at Q 2 ~ 5.7 and 7.0 
(GeV/c) 2 for centre-of-mass energies from threshold to 1.8 GeV, encompassing the Sn(1535) reso- 
nance, which dominates the channel. This is the highest momentum transfer measurement of this 
exclusive process to date. The helicity-conserving transition amplitude A1/2, for the production of 
the £11(1535) resonance, is extracted from the data. Within the limited Q 2 now measured, this 
quantity appears to begin scaling as Q~ 3 — a predicted, but not definitive, signal of the dominance 
of perturbative QCD, at Q 2 ~ 5 (GeV/c) 2 . 
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I. INTRODUCTION 
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The goal of strong interaction physics is to understand 
hadrons in terms of their fundamental constituents, the 
quarks and gluons. Although these constituents are de- 
scribed by Quantum Chromodynamics (QCD), and per- 
turbative methods work well where applicable, mostly 



the complexity of the theory precludes a description of 
hadrons in terms of QCD. Various techniques are used to 
make progress, such as numerical simulation of QCD and 
hadron models with effective QCD degrees of freedom. 
In this sense, the role of experiment is to make measure- 
ments which test the predictions of QCD-inspired quark 
models. Most models can describe the static nucleon 
properties and the baryon spectrum, and so other mea- 
surements, such as electromagnetic transition form fac- 
tors and strong decay amplitudes, are required. 

A baryon's quark substructure can be excited into 
a resonance — an excited state of the quarks with well- 
defined baryon quantum numbers. The transition form 
factor is the coupling (amplitude for the transition) from 
one baryon state to another, as a function of the squared 
invariant momentum transferred to the baryon Q 2 . The 
measurements of couplings between baryon states and 
the dependence of these on Q 2 , can be used as stringent 
tests of quark models. These couplings can be expressed 
in terms of the transition matrix elements between states 
of definite helicity. 

The difficulty in measuring baryon transition form fac- 
tors lies in isolating any of the multitude of wide and over- 
lapping resonant states. The 5n(1535) is a baryon res- 
onance that can be accessed relatively easily. Although 
there are many overlapping states in its mass region, it 
is very strongly excited over the accessible Q 2 range and 
is the only resonance with a large branching fraction to 
77 mesons ill], causing it to dominate the p(e, e'p)rj chan- 
nel. This dominance is partly due to isospin conservation, 
since the proton has isospin=i and the 77 has isospin=0, 
only the N*(I — |) resonances can decay to a proton- 77 
final state — N*(I = §) resonances are forbidden. 

As well as being accessible, the Sn is an interesting 
resonant state. It is the negative parity partner of the 
nucleon, they are both spin-half and isospin-half parti- 
cles. The transition form factor for the production of 
the Sn falls more slowly with Q 2 than the dipole form 
factor G D = (1 + Q 2 /0.71)- 2 , at least up to Q 2 = 3.6 
GeV 2 0, and more slowly than the form factor for typ- 
ical baryons. An example is the L>i3(1520) Q, which is 
from the same SU(6) (££)0(3) multiplet and mass region 
as the 5n(1535). The S , n(1535) branching fraction to 
prj, at b v ~ 50%, is anomalously high when compared to 
that of the other N* resonances, a phenomenon which is 
not well understood. 

It is expected from helicity conservation in perturba- 
tive QCD (pQCD) that at sufficiently high Q 2 the photo- 
coupling amplitude A 1 / 2 will begin to scale as l/Q 3 Q, 
or equivalently the quantity Q 3 Ai/ 2 will flatten. The ob- 
servation of such scaling is thus a possible signal of the 
transition to the dominance of hard processes. This mo- 
tivates the present experiment which studied exclusive 77 
production, allowing access to the amplitude A\j 2 for the 
Sii resonance, at the highest ever Q 2 yet measured. 

The first measurement of 77 production at substantial 
Q 2 was published by Brasse et al. [5| in 1984 based on 
work at DESY that went to Q 2 = 2.0 and 3.0 GeV 2 . 



This was the first indication that the Sii(1535) falls far 
slower with Q 2 than the £>i3(1520), and hence dominates 
the channel at high Q 2 . 

In 1999, Armstrong et al. [2| published data obtained 
in Hall C at Jefferson Lab at Q 2 = 2.4 and 3.6 GeV 2 , the 
highest until this work. The cross-section was found to be 
about 30% lower than the DESY data and the full width 
of the 511(1535) about twice as wide. By comparing with 
inclusive data, a lower bound was put on the branching 
fraction Sn — > rjp of b v > 0.45. 

A recent paper by the CLAS collaboration from Hall 
B at Jefferson Lab jg], published data for this process, 
at centre-of-mass (cm.) energy W — 1.5 — 2.3 GeV and 
Q 2 = 0.13 — 3.3 GeV 2 . The photocoupling amplitude 
A1/2 of the proton to £11 (1535) transition was extracted, 
and the anisotropics in the differential cross-section were 
more precisely determined. The results for the magni- 
tude and width of the Sn resonance favoured the Arm- 
strong data over the older Brasse result. Evidence was 
shown for a significant contribution to 77 electroproduc- 
tion due to a P-wave resonance with a mass around 1.7 
GeV. 

This paper describes an experiment where electrons 
were scattered off free protons at high momentum trans- 
fer and both electron and proton were detected in co- 
incidence. In Sec. |TT] the kinematics of the reaction are 
discussed along with the formalities of the cross section 
and helicity amplitude for 511(1535) production. Sec. IIIII 
describes the apparatus and methods used to acquire the 
data. Sec. IIVI then goes on to present the processing of 
the data including corrections, calibration, cuts, Monte 
Carlo simulation, backgrounds and ultimately the cross 
section extraction and error analysis. In Sec.|V]the 77 pro- 
duction differential cross section is plotted and fit with 
an angular dependence. A Breit-Wigner form is fitted to 
the data and the Sn helicity amplitude and resonance 
parameters are extracted. A brief summary is given in 
Sec. EH The appendix tabulates the extracted 77 produc- 
tion differential cross section. 



II. FORMALISM 

A. Kinematics 

Figure Q] shows the one-photon exchange (Born) dia- 
gram for the resonance electroproduction process. The 
incident electron ki scatters off the stationary proton 
Pi with mass m p . We detect the scattered electron kf 
and proton pf and reconstruct the undetected particle p x 
using the missing mass technique, evaluated from four- 
momentum conservation 



(pi + h -p f - kf) 2 . 



(1) 



Using the symbols from the diagram and neglecting the 
electron mass, the positive square of the four-momentum 



h = {E, k t ) 



Pi = (m p ,6) 



k, = (E',k f ) 




Px = (E x ,p x ) 



Pf = ( E 'p>Pf) 



B. Cross Section 

The five-fold differential cross-section for the reaction 
may be expressed as the product of the transverse virtual 
photon flux Tt and the centre-of-mass cross-section for 
the electroproduction of the pr\ pair 

^ =w ^ iw :* (2) 

where the flux of transverse virtual photons in the Hand 
convention Ml is 



FIG. 1: The one-photon exchange diagram of the resonance 
electroproduction process, where, for example, fe, is the four- 
momentum vector of the incoming electron composed of en- 
ergy E and momentum fc;. 



T T (W,Q 2 ) 



a W K 1 
4tt 2 m p E 2 Q2 1 - e ' 



(3) 



the longitudinal polarisation of the virtual photon is 
given by 



transferred from the lepton to hadron system is Q 2 = 
—q 2 = <lEE'sm 2 (6 e /2). The mass of the resonant state 
is W 2 = (q + p.) 2 = q 2 + m 2 p + 2m p v. 

Figure [5] shows the scattering and reaction plane coor- 
dinate systems: 9 e is the scattering angle of the electron; 
9 pq the angle between the outgoing proton and the mo- 
mentum vector of the virtual photon, q; the polar and 
azimuthal angles of the missing momentum are 9* and (f> x 
respectively, defined with respect to q and the electron 
scattering plane. A super-scripted * denotes measure- 
ment in the pr\ centre-of-momentum frame. 




l + 2J§!rtan 2 (6> e /2) 



(4) 



and the energy required by a real photon to excite a 
proton to a resonance of mass W is 



K 



W 2 -mi 
2m T) 



(5) 



FIG. 2: The scattering and reaction plane coordinate systems. 



The unpolarised virtual photon cross-section is writ- 
ten in terms of the transverse polarised virtual photon 
dax/dfl* longitudinal polarised virtual photon da^/dfl* 
and interference contributions, da^x/dQ* and daxT/dQ* 
in Eq. ([6]) . Each of these four individual components are 
expressed in terms of multipoles p, Q, where Ei±, Mi± 
and Si± are the electric, magnetic and scalar multipoles 
respectively; I is the orbital angular momentum, and ± 
indicates the total angular momentum via j = I ± i . If 
only terms with I < 2 and either of the dominant isotropic 
multipoles, E Q+ or So+, are retained, then Eqs. © are 
obtained [10| . 

The virtual photon cross-section is parametrised in 
terms of its angular dependence as Eq. |(5J|. The pa- 
rameters A — F are then given in terms of the truncated 
multipolc expansion by the Eqs. ©. 
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— %- = A + B cosfl* + C cos 2 9* + D sm6*cos(f>* + E cos0*sin(9*cos<A* + F sin 2 6'*cos2fi!>* 
dQ* 
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C. Helicity Amplitude 

The helicity amplitude is the matrix element that con- 
nects states of definite (the same or different) helicity. As 
such, it is a convenient measure of the coupling strength 
between states and can be used to fundamentally test 
quark models. The amplitudes are labelled by the virtual 
photon polarisation (either transverse A or longitudinal 
S) and the total 7JV helicity (^ or I). Spin-i resonances 
are therefore described only by A\i 2 and Si/2- 

The helicity amplitude A 1 / 2 , for the process j v p — ► 
S'ii(1535), can be obtained from the contribution of the 
Sii(1535) to the Eq_i_ multipole at the resonant mass 



W = W R , using [11|, [12 1 



.4 



1/2 



2tt 



|p;i«w-flW fl r fl 



m p K m p brj 



1/2 



|^)+(Wk)|. (10) 



This requires, not only isolating the 511(1535) from 
the other resonances and the non-resonant background, 
but further isolating the Eo + multipole from the other 
multipoles. In this case, for j v p — > r\-p at the 5n(1535) 
resonance mass, such an isolation is almost implicit in 
the measurement due to the dominance of the 5n(1535). 
Being an S'-wave resonance, implies a dominance of the 
isotropic multipoles — which has previously been seen in 
the data [2, ua, la, LL3|. So too, among the isotropic con- 
tributions it appears that the transverse multipole Eq + , 
dwarfs the longitudinal part Sq + 0, [g, LL4J • 

Doing a longitudinal/transverse (LT) separation re- 
quires measuring the cross-section for at least two values 



of e at the same Q 2 , which was not done in this exper- 
iment. Such separations performed in the late 1970's 
are consistent with no longitudinal component. Where 
R = olJot is the ratio of longitudinal to transverse 
cross sections, Ref [3 found R = 0.23 ± 0.15 at Q 2 
= 0.4 GeV 2 and Refj| found R = 0.22 ± 0.23 at Q 2 
= 0.6 GeV 2 and R = -0.16 ± 0.16 at Q 2 = 1 GeV 2 . 
Quark models [15| show this ratio decreasing with Q 2 . 
In angular fit to their recent data, Denizli [6| shows that 
the parameters of the P 1 (cos0*)cos<fi components fluctu- 
ate around zero and are consistent with zero within ex- 
perimental uncertainty. These parameters measure the 
daLT/d^l component of da/dfl, suggesting that the lon- 
gitudinal component is small — but since daur/dQ* is a 
sum of terms with possibly different signs, it is possible 
that <So+ is in fact comparable to Eq+. 

In this paper it is assumed that the longitudinal ampli- 
tudes are not significant for this reaction. The validity of 
this will become clear in the future when LT separations 
are done at high Q 2 . 

The cross-section can thus be written as depending 
only on the dominant Eq + multipole in the simple form 



da 

dn 



>i 



\p;\w 

m p K 



Eq^ 



The combination of Eq. (fTTj) and Eq. (fT0|) yields 



^i/ 2 (Q 2 ) = J5t^(Q 2 ), 



2m p &,j 



(11) 



(12) 



the helicity amplitude as a function of gr = cr(Wjj) [the 



total cross-section of the Sii(1535) resonance, measured 
at the resonance mass Wr.] 

The Eq+ multipole can be more reliably extracted from 
a fit to the angular dependence. Parameters A and C in 
Eqs. ^ share some common terms, and a simple cancel- 
lation yields Eq. (|13|) — although in the absence of an LT 
separation, it still must be assumed that Sq-\- is negligible. 



A + \C 



E, 



o+l 



\p;\w 

m p K 

WW 2 

E 0+ \ 



Q 2 



US', 



o+l 



m p K 



(13) 



Where possible in this work, the E 0+ multipole is ex- 
tracted using both methods, but for consistency with pre- 
vious analyses the final result is quoted from the method 
assuming isotropy. 



III. THE EXPERIMENT 



HallC 




FIG. 3: A plan view of Hall C showing the beamline, target 
and the SOS and HMS spectrometers which detected electrons 
and protons respectively. Figure from Ref. Q|. 



The experiment, measuring the unpolarised differential 
cross-section for the process p(e, e'p)ri, was performed in 
Hall C (Fig. E]) of the Thomas Jefferson National Ac- 
celerator Facility, during May and June of 2003. The 
Short Orbit Spectrometer (SOS) [lq |. a resistive QDD 
(quadrupole, dispersive dipole, anti-dispersive dipolc) 
spectrometer, was used to detect scattered electrons. The 
High Momentum Spectrometer (HMS) [17fl . with a su- 
perconducting QQQD configuration, detected the recoil 



protons. The r\ particles were identified using the missing 
mass method. 

Both spectrometers have a similar detector ensem- 
ble, including drift chambers for determining the track, 
scintillator arrays for triggering, an electromagnetic 
calorimeter for particle identification (P1D) and a thresh- 
old gas Cerenkov also for PID and tuned to differentiate 
between pions and electrons in the SOS. Figure [U show- 
ing the detector components, is representative of cither 
detector stack. 




Pb-glass 
Calorimeter 



DC1 



FIG. 4: A side view of the HMS detector stack, which is 
also representative of the SOS. The detected particles travel 
from left to right, encountering first the two drift chambers 
(DC) then the first two arrays of scintillators (SI) oriented in 
the X and Y directions, then the gas Cerenkov detector, the 
third and fourth scintillator arrays and finally the calorimeter. 
Figure from Ref. [^ 



The Jefferson Laboratory's superconducting radiofrc- 
quency Continuous Electron Beam Accelerator Facility 
(CEBAF) provides multi-GeV continuous-wave beams 
for experiments at the nuclear and particle physics in- 
terface [18| . The accelerator consists of two anti-parallel 
linacs linked by nine recirculation beam lines in the shape 
of a racetrack, for up to five passes. Beam energies up to 
nearly 6 GeV at 100 \xk and > 75% polarization are pos- 
sible. For this experiment, the incident electrons had the 
maximum available energy: E e — 5.500 GeV for most of 
the experiment and E e = 5.491 GeV for an 11 day period 
near the beginning. 

The target was liquid hydrogen maintained at a tem- 
perature of 19 K. The beam passes through 3.941 cm 
of liquid and through 0.12 mm of aluminium target cell 
walls on entrance and exit. The beam was rastered within 
a square of ± 1 mm to minimise density changes due to 
target boiling. A dummy target consisting of two alu- 
minium plates was used to simulate reactions within the 
target walls. 

The trigger for the experiment was a coincidence be- 
tween pre-triggers (or singles triggers) from both of the 



spectrometers. Both of the spectrometer pre-triggers 
were the requirement of a signal in three out of the four 
scintillator planes (SCIN) . In addition to the coincidence 
trigger, data were taken for singles triggers from both of 
the two spectrometers. This was pre-scaled according to 
the rate so as not to interfere with the coincidence trigger. 
This singles data allowed the monitoring of the luminos- 
ity and the electron detection efficiency. The elastic scat- 
tering events within the SOS were used to monitor the 
beam energy and the performance of the SOS magnets. 

Blok et al. [l9| is descriptive of the accelerator, beam 
monitoring equipment and current monitors, target ras- 
tcring system, beam energy measurement and cryogenic 
target. More detailed discussions are made of the two 
spectrometers, their detector packages, the trigger logic 
and data acquisition. Further references are provided for 
all covered topics, the interested reader is advised to con- 
sult that work. 

The electron spectrometer was fixed in angle and mo- 
mentum, thereby defining a central three-momentum 
transfer vector q for the virtual photon which mediates 
the reaction. Around this q vector is a cone of reac- 
tion products including the protons from the resonance 
decay of interest in this measurement. The "kinematic 
focusing" caused by the high momentum transfer of the 
reaction makes it possible to capture a large fraction of 
centre of mass decay solid angle in a spectrometer, as it 
comes out as a "narrow" cone in the lab. The proton 
spectrometer was stepped in overlapping angle and mo- 
mentum steps to capture as much of this decay cone as 
possible. 

The exact choice of kinematics was based on a compro- 
mise between maximising the Q 2 for the available beam 
energy and detecting the full centre-of-mass decay cone 
for the p(e, e'p)ir° reaction to the highest possible W. 
This reaction, which was measured concurrently is re- 
ported on by Villano [20 ]. The maximum central mo- 
mentum of the SOS, 1.74 GeV, required increasing #sos 
to increase the Q 2 , while the minimum HMS angle of 10.5 
degrees required decreasing Qsos to extend the full angu- 
lar coverage to higher W. At Osos = 47.5 degrees and the 
maximum SOS momentum, it was found that the kine- 
matic region from pion threshold to above the Su mass 
fell nicely within the best resolution region of the SOS 
spectrometer and full cos9* coverage was possible for the 
pn up to W = 1.4 GeV and pr] up to W = 1.6 GeV, 
as shown in Fig. [5] These SOS central parameters cor- 
respond to a virtual photon with momentum 4.51 GeV 
and angle 16.5 degrees and Q 2 ~ 5.8 (GeV/c) 2 at the Su 
resonance mass. 

In addition to data taken with these kinematics, it was 
decided to take a smaller set of data at even higher Q 2 , 
although the angular coverage would be incomplete. In 
this configuration the SOS was set with central momen- 
tum of 1.04 GeV and angle of 70 degrees, which gives 
a central virtual photon with \q\ = 5.24 GeV and angle 
10.8 degrees and Q 2 ~ 7.0 (GeV/c) 2 at the Su resonance 
mass. For the purposes of this paper, the first data set 
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FIG. 5: The W acceptance of the detector pair for the 
e(p, e'p)r\ reaction, in the lower-Q 2 configuration, as a func- 
tion of the laboratory scattering angle and momentum of the 
proton. The contours are constant W of the hadronic sys- 
tem, for an electron at 47.5 degrees and momentum of 1.74 
GeV/c, for the full range of 8* and cf> — and 180 degrees. 
The solid central contour is W = 1.5 GeV, from which they 
increase in steps of 100 MeV to the outermost at W — 1.9 
GeV. In practice, the angle and momentum bite of the SOS 
causes the contours to be much broader. Each black box is 
the acceptance of a particular HMS setting (^Table|I|. The 
alternate settings are offset by 1.5° and are a 4.7% increase 
in momentum, so that they are approximately centred on the 
points where the boxes join. 



will be called the 'lower-Q 2 ' configuration and the second 
data set, the 'higher-Q 2 '. The kinematic settings for the 
experiment are summarised in Table [T] 



Electron Arm 




Proton Arm 






Psos 
GeV 


#sos 
degrees 


PHMS 

GeV 




#HMS 

degrees 








4.70 








18.0, 


15.0 






4.50 f 




19.5, 


16.5, 


13.5, 


11.2 






3.90 




21.0, 


18.0, 


15.0, 


12.0 






3.73 f 


22.5, 


19.5, 


16.5, 


13.5, 


11.2 


1.74 


47.5 


3.24 


24.0, 


21.0, 


18.0, 


15.0, 


12.0 






3.10 f 


22.5, 


19.5, 


16.5, 


13.5, 


11.2 






2.69 


24.0, 


21.0, 


18.0, 


15.0, 


12.0 






2.57 f 


22.5, 


19.5, 


16.5, 


13.5, 


11.2 






2.23 




21.0, 


18.0, 


15.0, 


12.0 






2.13 f 




22.5, 


19.5, 


16.5, 


13.5 






4.70 










11.2 






4.50 










14.2 


1.04 


70.0 


3.90 
3.73 
3.24 








14.2, 


11.2 
11.2 
11.2 



TABLE I: The kinematic settings of the two spectrometers. 



Data were taken at a mean beam current of 92 /zA. 
The lower-Q 2 configuration was run for 6 weeks, totaling 
127 C of electrons through the target from which about 
50,000 r\ particles were identified from proton-electron 
coincidences, by missing mass reconstruction. Due to 
improved accelerator operation, the one week of running 
for the higher-Q 2 setting received 29 C of charge, but 
only about 2,000 r\ particles were reconstructed. 



IV. DATA ANALYSIS 

The raw data as recorded by the electronics were re- 
played offline to produce PAW or ROOT ntuples of cali- 
brated physics quantities. Corrections were made to the 
data for inefficiencies, dead times and accidental coinci- 
dences. The detector response was simulated using the 
Monte Carlo technique (including multiple scattering in 
the detector and nuclear reactions in the target walls) 
with one input model cross-section for the r\ production 
signal and another model for the multipion background 
processes, described in detail in Sec. IIVCI Using an it- 
erative procedure, a linear combination of the signal and 
background simulations was fitted to the data and the 
result used to refine the simulation input model, until 
the simulation in each bin matched the data with a mul- 
tiplicative factor of close to unity. 



Effect 


lower 


-Q" higher-Q" 


Proton absorption 




+4 ± 1% 


f Computer DT 


+ (1.0 - 


19.1)% +(1.8 - 10.9)% 


f HMS tracking 


+ (2.3 - 


14.3)% +(3.3 - 7.4)% 


f SOS tracking 


+ (0.3 - 


0.9)% +(0.2 - 0.8)% 


tElectronics DT 


+ (0.0 - 


2.4)% +(0.0 - 0.6)% 


* Random coincidence 


-(0.0 - 


7.6)% -(0.0 - 1.2)% 



TABLE II: Corrections applied to the data. For corrections 
applied ^run-by-run or * bin- by-bin, the range of the size is 
indicated in parentheses. 



and coincident l H(e,e'p) events. From these kinemati- 
cally overdetermined reactions, it was possible to check 
the momentum p and angles 9 and cf> (in-plane and out- 
of-plane relative to the spectrometer central axis respec- 
tively) for both spectrometers, and the beam energy E. 
A fit was done to determine what offsets to these quanti- 
ties most accurately produced the required values for the 
invariant hadronic mass, and missing mass and energy 
for the elastic scattering. In the case of SOS momen- 
tum (equivalently the dipole field), there is a saturation 
as the current is increased due to the resistive nature of 
the magnets. A field dependent correction was thus de- 
termined. These offsets, summarised in Table IIII1 were 
used in the replay of the present data. 



A. Raw Data to Physical Quantities 

The raw data from each trigger was stored onto tape. 
These data were "replayed" offline a number of times dur- 
ing the analysis, using the Hall C data reduction code, as 
the calibration of the detectors was improved. For each 
event, a list of calibrated event properties including posi- 
tion and angles of the track, timing and energy deposition 
information were determined. So too were quantities for 
the scattering including the centre-of-mass angles, invari- 
ant hadronic mass and the missing mass. For each run 
an ntuple of these event parameters was produced along 
with a file containing scaler information and calculated 
efficiencies and dead times for that run. 

The data were corrected on a run-by-run basis for 
these inefficiencies and dead times during the filling 
procedure — each event passing the cuts was filled into 
the histogram weighted by a run dependent correction 
factor. This included track reconstruction inefficiencies 
in the HMS and SOS spectrometers and computer and 
electronic dead times. A summary of all the corrections 
applied to the data is given in Table [TTJ 

The pion form factor (F^) experiment [19|,[2l| was con- 
ducted in the same suite of experiments as the current 
experiment and this work makes reference to some analy- 
ses reported there. A detailed description of the fitting of 
the reconstruction matrix elements for the spectrometers 
is included there. A number of offsets and corrections 
were determined by analysing singles elastic scattering 



Quantity 



HMS 



SOS 



p ( lower- Q 2 ) 
p (higher-Q 2 ) 



0.0 ± 0.5 mrad 
+1.1 ± 0.5 mrad 
-0.13 ± 0.05% 



0.0 ± 0.5 mrad 
+3.2 ± 0.5 mrad 
-1.36 ± 0.05% 
0.00 ± 0.05% 



0.00 ± 0.05% 



TABLE III: Nominal 2003 spectrometer offsets [IE |2l| ap- 
plied to the data during the replay phase. 



1, Trigger Efficiency 

The HMS trigger was a three out of four coincidence 
between the four scintillator planes. A trigger inefficiency 
for proton detection in the HMS is produced by protons 
which are not detected in their interaction with the scin- 
tillator, and by protons that do not make it through all 
the scintillators due to absorption. 

A previous study of general HMS trigger efficiency (19 . 
|21| showed a strong dependence on relative particle mo- 
mentum i5hms- The momentum in the spectrometers 
is measured relative to the central momentum p so t, so 
that particles with the same S = (p — p sc t)/Psot are dis- 
persed by the same amount. The trigger efficiency was 
mostly very high at 0.995 but dropped rapidly for mo- 
menta lower than S ~ —6%. The data was analysed with 



a cut of 6 > —6 resulting in an average increase in ex- 
tracted cross-section of 1.4%. No correction for this effect 
was made, but this figure was used as an estimate of the 
error due to the trigger efficiency. 

The trigger requires hits in scintillator planes SI 
and S2, so another source of inefficiency is absorption, 
through nuclear reactions, of the proton in target or de- 
tector materials before the S2 plane. The total pp colli- 
sion cross-section, a pp , varies slightly from 47 to 42 mb 
for proton lab momenta between 2 to 5 GeV/c [lj which 
is the momentum range of this experiment. Therefore for 
this experiment, the trigger efficiency due to absorption 
is relatively independent of kinematic setting. 

The primary sources of interacting material are the SI 
scintillator planes which had a thickness of 1 cm each and 
the Aluminium windows of the gas Cerenkov and aerogel 
detectors which had a total thickness of 0.51 cm. The 
proton- nuclear cross-section was estimated as A or a pp . 
Combining interactions in all material, the trigger effi- 
ciency due to proton absorption is estimated to be 0.95. 

To calculate the correction used in the experiment for 
the trigger efficiency due to proton absorption, a study 
of ep elastic events was done. The SOS was set for elec- 
trons at central angle = 50° and central momentum of 
1.74 GeV/c and the HMS was set for protons at central 
angle of 18° and central momentum of 4.34 GeV/c at a 
beam energy of 5.247 GeV. For a point target, the SOS 
has an out-of-plane angular acceptance of ±37 mr and 
an in-plane angular acceptance of ±57 mr (the in and 
out-of-plane angles are relative to the central axis of the 
spectrometer), while the HMS has an out-of-plane angu- 
lar acceptance of ±70 mr and an in-plane angular accep- 
tance of ±27 mr. The ratio of electron to proton momen- 
tum is 0.4, so for the maximum SOS out-of-plane angle, 
the corresponding HMS out-of-plane angle is 15 mr. The 
maximum SOS in-plane angle gives a corresponding HMS 
in-plane angle of 23 mr. The data acquisition is set-up 
to accept singles triggers from the SOS and HMS indi- 
vidually in addition to coincidence triggers between the 
HMS and SOS. 

In offline analysis, the cuts described in Sec. IIVB II 
were used to identify electrons in the SOS. A good elastic 
event in the SOS was identified by a cut of 0.9 < W < 
1.0 GeV and a cut on the SOS in-plane angle of ±50 mr 
which ensured that the proton would be within the HMS 
angular acceptance. The proton was selected by the time- 
of-flight between electron and the particle detected in the 
HMS. The raw number of coincident ep events was 2009 
and the number of single events was 205. 

Some data were taken with an aluminium "dummy" 
target, which is intended to model an empty target cell, 
but is 7.78 times thicker in order to increase the count 
rate. Analysis of this data determined that the target 
endcaps would contribute 12 ± 3 events to the raw co- 
incidence events and 123.0 ± 12 events to the raw single 
events. Therefore the proton trigger efficiency due to 
absorption is 0.96 ± 0.01, in good agreement with the 
prediction. A 4% correction was applied to the data for 



this effect. 



2. Calibration of Simulation Resolution 



As described in Section [TVB31 the cross section is ob- 
tained by integrating over the r\ missing mass peak. The 
sensitivity to the exact shape of the peak is thus small 
and is minimised by matching the simulation resolution 
to the data resolution as far as possible. This was done 
for the elastic peak, following which the simulation of the 
eta peak was in agreement with the data width without 
further change. 

The elastic scattering of electrons into the SOS and 
protons into the HMS were compared to SIMC Monte 
Carlo simulations of the same. The invariant mass deter- 
mined from elastic scattering must be the proton mass, 
but is broadened due to resolution effects and radiative 
tails — which are included in the simulation. In both de- 
tectors it was found that the width of this peak predicted 
by the simulation was narrower than for the data. These 
resolution differences were taken into account by increas- 
ing the drift chamber resolutions in each spectrometer, 
from the nominal value of 300 /im. The resolution was 
varied, and the simulation repeated, until a Gaussian fit- 
ted to the simulated spectrum had the same width as 
a Gaussian (and polynomial background) fitted to the 
data. It was found that the HMS needed a drift cham- 
ber resolution of 570 /im to match the data width of 
17.2 MeV, while in the SOS, the 30.1 MeV width was 
achieved with a 350 /im resolution. 

This method had the effect of degrading the optics of 
the simulation slightly to match the experimental trans- 
port matrix elements of the data, in a logical yet simple 
manner. During the systematic error analysis process, 
described in Section TlVFl the drift chamber resolutions 
were varied by 10% and the integration limits were var- 
ied by 0.1 GeV 2 , characterising the sensitivity of the ex- 
tracted differential cross-sections and amplitudes. 



3. Collimator Punch Through 

A source of background is due to particles that interact 
with the edges of the HMS collimator aperture, located 
just before the first quadrupole magnet, whose kinemat- 
ics are thus changed. The collimator is made from 6.35 
cm thick HEAVYMET (machinable Tungsten with 10% 
CuNi; density=17 g/cm 3 .) For practical purposes elec- 
trons are stopped by the SOS collimator, but protons 
have the possibility of "punching" through the collima- 
tor, undergoing multiple scattering and energy loss in the 
material, and still making it through the spectrometer to 
the detectors. This process is modelled in the simulation 
of the experiment and additionally a loose cut, 3 cen- 
timeters outside the collimator edge, is used to eliminate 
unphysical reconstructions. 



4- Data Cuts 

The 'standard' cuts are listed in Table IIVI The cuts 
on relative electron momentum <5sos > and relative proton 
momentum <5hms, are made to ensure that only parti- 
cles within the well understood region of the spectrom- 
eter momentum acceptance are used. The momentum 
in the spectrometers is measured relative to the cen- 
tral momentum p so t, so that particles with the same 
S = (p — p S et)/Pset are dispersed by the same amount. 

Some parts of the SOS spectrometer acceptance, due 
to an ambiguity in the solution of the optics equations, 
do not reconstruct reliable tracks. The cuts on the SOS 
focal plane position in the magnet dispersion direction, 
^SOS,f.p., are to eliminate these regions. The particle 
identification cuts are described in Section IIVB1 



Quantity 



Variable Cut 



Electron momentum 


<feos 


< +20% 
> -15% 


Proton momentum 


5hms 


< +9% 
> -9% 


SOS focal plane position 


^SOS,f.p. 


> -20 cm 


dispersive direction 




< +22 cm 


' Coincidence time 


I^coin £cent | 


< 1.5 ns 


tSOS Cerenkov 


iVp. e . 


> 0.5 


^SOS calorimeter 


-C'liorm 


> 0.7 



TABLE IV: The set of 'standard' cuts applied to the data and 
to the simulations where applicable. ^The Particle Identifica- 
tion cuts are not applied to the simulation. 



5. Binning 

The data were binned in W, cos6>* <f>* and m 2 , where 
W is the invariant mass of the hadronic system, 9* is 
the polar angle between the direction of the 77 and the 
three-momentum transfer vector q in the centre-of-mass 
of the resonance, cf>* is the azimuthal angle of the r\ with 
respect to the electron scattering plane, and m 2 is the 
square of the missing mass for p(e, e'p)x. 

For the lower-Q 2 data, this was done in 12 cos#*-bins 
and 8 </>*-bins, to maximise the angular resolution for 
partial- wave analyses, necessitating m 2 -bins of 0.1 GeV 2 
and IF-bins of 30 MeV near the resonance and 40 MeV at 
higher W . The higher Q 2 data, with far fewer detected 
particles, was binned with VT-bins of 30 MeV, 6 cos6>*- 
bins, 5 0*-bins and m 2 -bins of 0.15 GeV 2 . 

Bins in (W, cos9* </>*) were retained for the analysis 
if they passed the following three criteria. Firstly, in the 
region of the r] missing mass peak, the simulation was re- 
quired to predict a signal to background ratio of at least 
0.25. Secondly, the simulation needed to have predicted 
a minimum average number of 77 events in the peak of 
1.5 per missing mass squared channel. This criterion was 



used instead of requiring a total number of predicted r\ 
particles because the resolution of the missing mass peak 
changes substantially with cos6>* The third criterion for 
acceptance was, following the subtraction of the all the 
backgrounds, the sum of the data in the region of the 
missing mass peak was required to have a statistical un- 
certainty of less than 50%. 



B. Particle Identification 

1. Electron Identification 

In the SOS spectrometer, the Cerenkov detector and 
the electromagnetic calorimeter were used to identify 
electrons and reject pions. The Cerenkov detector was 
filled with Freon-13 at 1 atmosphere, yielding a veloc- 
ity threshold of (3 t = 1/n = 0.9992. The highest mo- 
menta detected by the SOS in this experiment was about 
2.09 GeV/c, corresponding to (3 = 0.9978 for pions, which 
is below the threshold for detection while all electrons are 
well above the threshold. Some pions make small signals 
in the Cerenkov due to scintillation or "knock-on" elec- 
trons from atomic scattering. The detected signal was 
calibrated into units of the number of photo-electrons, 

iV p . e .. 

For each event, the signals from each of the 44 lead- 
glass blocks in the calorimeter were summed to obtain the 
total energy deposited, E to t- This energy was then nor- 
malised by the momentum of the particle as determined 
by the tracking, ptrack, to obtain E noTm = -Etot/ptrack- 
The 16 radiation lengths of lead-glass bring electrons to 
a stop, resulting in a peak at -E nor m ~ 1 due to elec- 
trons. The pions peak at about -Enorm ~ 0.25, but have 
a long tail to higher E noIm due to the charge exchange 
nuclear interaction ir~p — > ir°nx, and subsequent decay 

7T° — * 77. 

Figure[S]shows the correlation between -E n orm and iV pe . 
for the lower-Q 2 data. The electrons are clearly well sep- 
arated from the pions by these two detectors. In the 
analysis, electrons are identified using two simple cuts, 
iVp. c . > 0.5 and E noim > 0.7, shown in the figure. 



Proton Identification and Accidental Coincidence 
Subtraction 



Protons were separated from pions using time of flight 
considerations. The raw difference in arrival times, idiff , 
between the electron in the SOS and the positive particle 
in the HMS, were corrected event-by-event for differences 
in path length of both particles through the detectors 
and the variation in velocity /3 of the positive particle 
(all electrons having essentially the same velocity.) This 
corrected coincidence time, i C oin, is plotted in Fig. [7] and 
shows peaks due to protons and n + particles and a back- 
ground of accidental (or random) coincidences. 
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FIG. 6: The correlation between E n 



and iV D .o. for all 



the lower-Q data. The particle ID cuts to select electrons, 



AS, 



> 0.5 and E noTrn > 0.7, are visible as dashed lines in 



the figure. All other cuts listed in Table HV] have already been 
applied to the data. 



FIG. 7: Coincidence time spectrum for the lower-Q 2 data with 
all of the 'standard' cuts except the coincidence time cut. The 
dark grey shaded region represents the 3 ns wide proton cut. 
The 2 ns beam structure is clear in the accidental background. 
Data from the light grey shaded regions was used to estimate 
the amount of accidentals in each (W, cosO^, 0JJ and m 2 ,) bin 
under the proton peak. 



The path taken is determined by the tracking algo- 
rithm from drift chamber hit positions while the velocity 
(3 = (p 2 /(m 2 + p 2 )) 1 / 2 is calculated from the measured 
momentum p assuming the proton mass m p . For pro- 
tons the corrected coincidence time depends only on the 
actual difference in starting times of the particles in the 
target, causing a peak of real coincidences, which has 
been shifted to zero in the figure. Particles with a dif- 
ferent mass, such as pions, have their coincidence time 
peak shifted relative to the protons since for the same 
momentum, they have a different velocity. The 7r + peak 
is broader than the proton peak because t coln is calcu- 
lated to remove the momentum dependence of the pro- 
tons but the pion locus remains momentum dependent. 
A much smaller number of kaons are detected and form 
a locus between the pions and protons, but remain dis- 
tinctly separable. It was then possible to select the pro- 
ton events and reject the pion and kaon events and most 
of the accidental coincidences using one simple cut. 

Accidental coincidences occur when both detectors are 
triggered within the 100 ns coincidence time window, 
but the detected particles originate in different scatter- 
ing events. In the coincidence time spectrum of Fig. 
the accidentals are the continuous background under the 
two main peaks. The 2 ns beam structure can clearly be 
seen in the spectrum. A 3 ns particle identification win- 
dow was used to select protons, but within this cut there 
is still some background due to accidental coincidences 
which must be subtracted. 

For each bin in 4 dimensions (W, cos6>*, <j>* and m 2 ), 
the number of accidental coincidences inside the proton 
cut was estimated by determining the average number 
of accidentals in the "wings" of the spectrum, —45 ns < 
tcoin < ~5 n s and 15 ns < £ co i„ < 50 ns, away from loci 
for actual coincidences. This value was then normalised 



for the width of the proton cut and subtracted from the 
data. The accidental correction is small for our kinemat- 
ics and rates, the weighted mean correction was 1.5% 
and the largest correction in any (W , cos6*, </>*) bin was 
7.6%. 



3. rj Identification 

In the case of inelastic scattering, the detection of the 
scattered electron and recoil proton is not an exclusive 
measurement — there will be at least one other emitted 
particle. If there is only one undetected particle it is pos- 
sible to fully reconstruct the kinematics of that particle. 
The data corresponding to such a channel, the p(e, e!p)r\ 
in this case, is isolated by constructing the square of the 
missing mass m 2 , as given in Eq. (Q]). 

Figure [5] shows the m 2 distribution for the lower-Q 2 
data, with the 7r°, rj and u> products are visible as peaks. 
The energy calibration of the data is done on the elastic 
peak. A simple Gaussian and polynomial background 
fitted to the eta peak is centered at 0.2993 GeV 2 , within 
0.2% of the nominal eta mass. The broadening is due to 
instrumental resolution. 

The actual extraction of the r\ particles is done by ap- 
plying a cut on m 2 around the 77 peak and subtracting 
the background. The resolution of this peak varies as a 
function of cos9* and therefore so does the cut, which is 
listed in Table El for the lower-Q 2 data. The higher-Q 2 
data has very little coverage above cos6** = at any W, 
and larger m 2 bins, so the cut was kept at a constant 
0.255 GeV 2 < m 2 < 0.36 GeV 2 . 

The continuous background, seen in Fig. [5J is due to 
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FIG. 8: Missing mass squared m 2 , from Eq. {T}, for all the 
lower-Q 2 data. 



events with more than one undetected particle. In this 
case, the missing mass does not correspond to any phys- 
ical mass because the magnitude of the missing momen- 
tum is smaller than the sum of the magnitudes of the 
individual momenta of the undetected particles. This 
effect, predominantly due to the production of multiple 
pions, is the principle background in this experiment, and 
is treated in Section IV C 3 ; 



"~~^^ cos#^ 


-0.917 


-0.750 


-0.583 


-0.417 


-0.250 


-0.083 


m 2 min (GeV 2 ) 
m 2 max (GeV 2 ) 


0.27 
0.34 


0.27 
0.35 


0.27 
0.35 


0.26 
0.36 


0.26 
0.36 


0.25 
0.37 


~---~^^ cos#^ 


0.083 


0.250 


0.417 


0.583 


0.750 


0.917 


m 2 min (GeV 2 ) 
m 2 max (GeV 2 ) 


0.25 
0.37 


0.25 
0.38 


0.25 
0.38 


0.25 
0.39 


0.25 
0.39 


0.25 
0.39 



TABLE V: The m x cuts used in each cos#* bin for the lower- 
Q 2 data. 



C. Monte Carlo Simulation of the Experiment 



11 



SIMC as a package consists of an event generator, 
which is able to produce events from a variety of physi- 
cal scattering processes common in Hall C or from phase 
space, and two 'single arm' spectrometer models, one 
for each detector, to track the particles and determine 
whether they are accepted by the detector. In each spec- 
trometer model, the particle is propagated from its initial 
position in its initial direction with transport maps pro- 
duced by COSY Infinity [23|], an arbitrary-order, beam 
dynamics simulation and analysis code, using the results 
of a field map of the magnetic elements. At points where 
there are apertures in the spectrometer such as collima- 
tors or the magnets themselves, the positions of the par- 
ticles are checked against these. For the magnets this is 
done at the entrance, exit and at the maximum beam en- 
velope within the object. Particles making it into the de- 
tector hut underwent multiple scattering and energy loss 
in the air and other materials. Particles that didn't con- 
form to the experimental trigger, such as passing through 
three of the scintillator hodoscopes, and for electrons the 
Cerenkov and calorimeter, were considered undetected. 
Detected events were reconstructed back to the target 
using the COSY optics matrix. 

SIMC was not used 'out of the box' for the present 
analysis, as it didn't have physics models for either the 
p(e,e'p)r? process or for multiple pion production. For 
r\ production, a simple model of the 5*11 resonance was 
added to SIMC, which was then run to simulate the signal 
part of the experiment. In the case of the multipions, an- 
other event generator was used and the resulting electron 
and proton pairs were propagated through the SIMC de- 
tector models to simulate their detection. Both the elas- 
tic ep and the ir° peaks have radiative tails that extend 
into the region m 2 > 0.1 GeV 2 . A Monte Carlo simula- 
tion determined that this makes a negligible contribution 
to the r\ peak compared to the multipion background. 

For both of the two Q 2 configurations, the data are 
taken in "settings" for which the HMS spectrometer an- 
gle and momentum is fixed. To limit file sizes and aid 
in online checking of the data, the data in each setting 
is taken in a number of "runs" . The simulation is per- 
formed on a run-by-run basis to match the data. The 
data and simulation are then binned into identical four- 
dimensional histograms. 



The Monte Carlo simulation of the experiment was 
done with SIMC [HI, the Jefferson Lab Hall C in- house 
detector simulation package. The simulation includes de- 
tailed models of both magnetic spectrometers and simu- 
lated the effects of radiative processes, multiple scatter- 
ing, and ionisation energy loss (due to material in the 
target and spectrometers). It was used to obtain the ex- 
perimental acceptance and radiative corrections for the 
resonance process under study, to simulate the multipion 
background to the resonance production and to study a 
number of other processes serving to verify our under- 
standing of the apparatus. 



1 . Model for p(e, e 'p)r) 

The model for r\ production used in the simulation 
and extraction of the cross-section is a single relativis- 
tic Breit-Wigner shape as a function of W multiplied by 
a exponential form factor depending on Q 2 . The form 
used for the Breit-Wigner resonance shape (from Christy 
and Bosted [24|) is given by 

lf_ if cm 

BW(1T) 



K(W)K cm (W) 
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T^totp7 



, (14) 



r[(w 2 -w^) 2 + (w R r tot ) 2 i 

where the equivalent photon energy in the lab frame is 



K{W) 



(W 2 - ml) 
2m„ 



the equivalent photon energy in the center of mass (CM) 
frame is 



K cm {W) 



{W 2 - ml) 



2W 



and Kr and K R m represent the same quantities evaluated 
at the mass of the Su resonance, Wr. T tot is the full 
decay width defined by 



r**^/^, 



(15) 



with (3j the branching fraction to the j decay mode and 
Tj the partial width for this decay mode. The partial 
widths are determined from the intrinsic widths T, using 



r\ = r 



p) 



pT\w r 



2L + 1 



(fi m \w R ) 2 +X* 



(16) 



where p^" 1 is the momentum, in the center of mass, of 
a meson produced by a system of invariant mass W and 
p°: m | yy R is the momentum of a meson from a system of the 
nominal resonance invariant mass Wr, L is the angular 
momentum of the resonance, and X = 0.165 GeV is an 
empirical damping parameter. The model as used in the 
simulation is then given by 



da 
dQ* 



1 

47r' 



-ae 



BW(W). 



(17) 



Although simplistic, the model describes the data well. 
The parameters a and b were obtained by fitting the form 
ae~ b ® to the cross-section at the Sii resonance mass, 
ctr, of data taken by Armstrong et al. Q and both of the 
present Q 2 data sets. The parameters Wr and Fr were 
refined using an iterative procedure in which the Breit- 
Wigner form was fitted to the angle-integrated lower-Q 2 
data, used to extract a new cross-section and then refit- 
ted. There was no explicit cos8* or (f>* dependence in the 
input model since the data showed very little anisotropy. 
The final model parameters are given in Table IVII 



2. Model for Multipion Production 



Parameter 


Value 


a 


9.02 nb 


b 


-0.479 (GeV/c)" 2 


Wr 


1525 MeV 


r„ 


133 MeV 


X 


0.165 GeV 



TABLE VI: The parameters of the Su resonance-dominated 
cross section model used for the final data extraction. 



the events should be generated. Depending on the re- 
action, the events are then sampled from interpolated 
data tables or according to a cross-section model — in 
contrast to SIMC behaviour which throws events uni- 
formly and weights them event-by-event. The generator 
itself extrapolates the cross-section from where data ex- 
ists to higher Q 2 using the square of the dipole form, 
(1 + Q 2 /0.71)~ 4 . The reactions included in our simula- 
tion of the multipion background are given by Eqs. (fl8|) 
and |H]). 



P 



P 



(model) (18) 



e+p^ e' + p + ir + ir~ it 

e+p^ e' +p + it + tc~it + tt~ (tables) (19) 

e+p^ e' + p + it + ir~ it + ir~ ir° 

The event generator was developed from an initial ver- 
sion for real photons [25j . In that version, for perfor- 
mance reasons, the cross-section is drawn from tabulated 
data — either measured or generated from models in un- 
measured regions. In the current version, the ptt + 7t~ 
exit channel, Eq. (|18|l . is now generated according to a 
phenomenological model [26j , with parameters that have 
been fit to recent CLAS data [27| which measured the 
process ep — > e'pit + ir~ for 1.4 < W < 2.1 GeV and 
0.5 < Q 2 < 1.5 GeV 2 /c 2 . The model is calculated for 
the three intermediate channels ir~A ++ , 7r + A° and pp. 
The amplitude is defined in the meson-baryon degrees 
of freedom, and is therefore not necessarily valid at this 
high momentum transfer Q 2 < 7 (GeV/c) 2 , where quark- 
gluon degrees of freedom may be the most appropriate. 
Radiative corrections are not implemented for the mul- 
tipion model. Despite these last two points, the results 
obtained are good enough to justify our implementation 
here. The properties of the generated pions aren't used, 
just the electron and proton pairs are propagated through 
SIMC. 



The multipion background was simulated using an 
event generator from the Jefferson Lab Hall B (CLAS 
detector) simulation package, which takes as input the 
Q 2 and W 2 ranges of the generation region and the re- 
actions, chosen from a list of possibilities, from which 



3. Multipion Background Subtraction 

As can be seen in Fig. [51 the peak at m 2 ~ 0.3 GeV 2 , 
corresponding to missing r\ particles, lies on a continuous 
background described in Sec. IIVB31 This background 
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was treated by simulating the m 2 spectra of the back- 
ground using SIMC with a model of the largest contribut- 
ing reactions, described in the previous section, and then 
subtracting the simulation from the data. 

The output of the simulation was a large set of multi- 
pion events that are accepted into our detectors. These 
events are then filled into histograms of the same struc- 
ture as those of the data, yielding our approximation to 
the shape of the multipion background, without an abso- 
lute normalisation. Since an absolute multipion cross- 
section is not being extracted, the shape is sufficient 
to subtract it from the data. It can be seen in Fig. [£l 
showing the data and associated simulations for one (W, 
cos6>*) bin, that excellent agreement is obtained. 

The simplest way to normalise the background to the 
data is with a two-parameter fit in each (W, cosO* <j>*) 
bin. The m 2 spectra of the multipion background simu- 
lation and the r\ production simulation would have been 
normalised to minimise the \ 2 difference between their 
sum and data m 2 spectrum. In practice, due to dimin- 
ishing acceptance, the out-of-plane <fr* bins demonstrate 
a phenomenon where the multipion background simula- 
tion and the r\ production simulation can have m 2 spec- 
tra similar enough to make a two-parameter fit unreli- 
able. This is typically the case for mid to large cos#* and 
worsens as W increases. An example of such a case is 
illustrated in Fig. [TO] 

For this reason, the fit was constrained to have the 
multipion normalisation parameter constant over <j)* as 
expected physically. For each and all of the (W, cos9*) 
bins, the fit had 9 parameters: one for the single mul- 
tipion normalisation over all the (f>* bins and one for rj 
production in each of the eight <f>* bins. The produc- 
tion of 7t° particles, seen as a peak at m 2 ~ 0.02 GeV 2 
in some panels of Figs. [TU] and IT2J produces a radiative 
tail which, in principle, extends under the rj peak. The 
size of this effect is smaller than the uncertainty in the 
multipion background, and so was neglected. 

This approach does a good job of reproducing the 
shape of the measured m 2 spectra. By eye, the sum of 
the normalised simulations seem to match the data well 
and in 94% of bins have a reduced \ 2 of less than 2. A few 
representative spectra showing the W and cos0* depen- 
dence of the m 2 distributions are shown in Figs.llllandfl2l 
respectively. The uncertainty in the normalised back- 
ground simulation was determined by adding the small 
Monte Carlo statistical uncertainty to the Minuit [28| fit 
uncertainty on the normalisation parameter in quadra- 
ture. 

It should be noted that some structure is seen within 
the normalisation parameters of the background model in 
W and cos#*, illustrated in Figs. IT51 and [T4l respectively. 
The extracted fit parameters seem to rise smoothly and 
approximately linearly with both increasing W and in- 
creasing cosO* . This is understandable since the multip- 
ion background model is produced from data with much 
lower Q 2 . Overall, the variation in the parameters is 
about a factor of 4. 
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FIG. 9: (Colour online) The <j> v dependence of missing mass 
squared distributions for W — 1.5 GeV and cosO^ = —0.916. 
The (green) points are the data while the solid fill is the sum of 
the simulations. The multipion background component is the 
light grey filled area and rj production simulation component 
has the darker grey fill. The dot-dashed lines shows the region 
within which the background fit is done while the dashed lines 
show the region within which the rj cross-section is extracted 



4- Target Window Background 

No explicit subtraction for scattering off the aluminium 
walls of the target was performed. The data taken with 
the dummy target in this experiment has too low statis- 
tics to be used for subtraction, and it was not taken at all 
of the experimental settings, but it is adequate for esti- 
mating the yield from the target walls and demonstrating 
the shape of the missing mass distribution. 

The size of the target wall effect is small and it has 
a very similar shape to the multipion background, so it 
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FIG. 10: (Colour online) The <j> v dependence ol missing mass 
squared distributions for W = 1.5 GeV and cosO^ = 0.416. 
Symbols as in Fig. \W\ Panels with fa = 1.178, 1.963, 4.320, 
and 5.105 are the out-ol-plane <f> bins where the simulations 
of the signal and background are sufficiently similar to make 
a two-parameter bin-by-bin fit unreliable. 



is therefore adequately accounted for in that background 
subtraction procedure. To first order, a nucleus is a bag 
of nucleons, and as such the multipion production from 
the aluminium target window has the same broad kine- 
matic distribution as from a free proton — the following 
analysis confirms this. 

The dummy target produced 430 coincidences from a 
beam charge of 1.97 C giving an average yield, integrated 
over all angles and W up to 1.7 GeV, of about 0.2 counts 
per mC. The hydrogen target's 64,000 multipion coinci- 
dences, estimated from the background subtraction pro- 
cedure, came at about 0.6 counts per mC, or three times 
as fast. Taking into account the differences in thick- 



ness between the dummy and the actual target walls, 
the multipion background is expected to have produced 
at least 20 times more background events than the tar- 
get walls. Figure [15] shows the similarity between missing 
mass spectra of the dummy data and the multipion back- 
ground simulation in three W bins. 



5. Radiative Corrections 

Radiative effects occur because photons are emitted 
in the interaction of the incoming and outgoing charged 
particles of the scattering. These real photons arc either 
produced within the field of the scattering nucleus itself, 
called internal radiation, or from the fields of other nu- 
clei in the propagation medium, called external radiation. 
This radiation causes there to be a difference between the 
actual momenta of the particles at the scattering vertex 
and the detected momenta, leading to measured values 
of W, Q 2 and the cm. angles cos9* and cf>, different 
from that of the actual scattering. In order to extract 
meaningful information from the detected particles, this 
radiative contribution must be corrected for. 

External radiation is small for the proton due to its 
high mass and can be handled essentially exactly for the 
electron, both pre and post-scattering. Dealing with in- 
ternal radiation requires a knowledge of the coupling of 
the photon to the electron, which is well known, and to 
the proton, which isn't known analytically since it de- 
pends on its QCD structure. It is then further compli- 
cated by interference of the amplitudes for radiation from 
each of the particles of the scattering. The radiative cor- 
rections for this experiment are done within SIMC, with 
the formalism of Ref. [291 ] , which is a general framework 
for applying radiative corrections in (e, e'p) coincidence 
reactions at GeV energies. This approach uses the angle 
peaking approximation and takes into account higher- 
order bremsstrahlung effects, multiple soft photon emis- 
sion and radiation from the scattered hadron. External 
radiation is also included in the model. 

The size of the radiative corrections implemented by 
SIMC is determined by running the full simulation with 
and without including radiative effects. In each bin, the 
ratio of the number of events predicted by these two sim- 
ulations, after the 'standard' cuts of Table IIVI and the 
missing mass cuts of Sec. IIV B 31 gives a number equiva- 
lent to the correction factor required to take account of 
the radiative effects. This radiative correction factor is 
listed for each bin in Tables IXIII and IXIIII along with 
the extracted cross-sections. Using these values and the 
size of the missing mass cuts given in Sec. IIV B 31 one 
can remove the effect of the radiative corrections on the 
cross-sections. 

The correction factor is plotted for the lower-Q 2 config- 
uration as a function of 4> for different W bins and three 
cos6>* ranges in Figs. [TUl [T71 and \TE[ The points are plot- 
ted for the kinematic bins where the data are sufficient 
to extract a cross section. Much of the large kinematic 
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FIG. 11: (Colour online) The W dependence of missing mass squared distributions for cos#* 
Symbols as in Fig. 1101 
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FIG. 12: (Colour online) The cos#* dependence of missing mass squared distributions for W = 1.5 GeV and <f) — 3.534 radians. 
Symbols as in Fig. [TO] The vertical dashed lines, which show the region within which the rj cross-section is extracted, vary with 
cos6** to accommodate the changing resolution. 
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FIG. 13: (Colour online) The W dependence of the normali- 
sation coefficient of the multipion background simulation. 



FIG. 14: (Colour online) The cos#* dependence of the nor- 
malisation coefficient of the multipion background simulation. 



dependence in these plots comes about due to the lim- 
ited acceptance, which decreases with increasing W and 
cos#*. 

This approach does neglect 2-photon radiation, which 
is expected to be about a factor of a = 1/137 smaller, and 
makes approximations. The uncertainty in the radiative 



corrections was estimated to be 2%. 
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FIG. 15: (Colour online) The (green) points are the m^ dis- 
tribution of the full set of data off the "dummy" target cell, 
and the grey filled histogram is the simulation of multipion 
background from Hydrogen. The simulation is arbitrarily nor- 
malised to match the data, with the same factor in all three 
panels. Note the similarity in shape. 
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FIG. 16: (Colour online) Radiative corrections for -1 < cos9 v 
< -|. Uncertainty is due to Monte Carlo statistics only. 



FIG. 17: (Colour online) Radiative corrections for -| < cosf?,, 
< | . Uncertainty is due to Monte Carlo statistics only. Sym- 
bols as in Fig. 1161 
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D. Extraction of the r\ Differential Cross-Section 

The actual r\ cross-section extraction is done by com- 
paring the data, having had the randoms and multip- 
ion background already subtracted, with a Monte Carlo 
simulation of the experiment, produced using SIMC and 
the T] production model described previously. The com- 
parison is done for each (W, cos9*, </>*) bin. The m x 
dependence of both the subtracted data and the simula- 
tion is integrated out between two tight limits in m x that 
contain the missing r\ particle peak 



Jhigh 

N l = J2 Nij i 

Jiow 



where i labels the (W,cos8* <f>*) bins and j labels 
the m x bins so that A™ is the content of a certain 
(W, cos8* 4>*m x ) bin. The number of simulation events, 



FIG. 18: (Colour online) Radiative corrections for | < cos8 v 
< 1. Uncertainty is due to Monte Carlo statistics only. Sym- 
bols as in Fig. 1161 



iV^j C , is obtained by multiplying the yield output of 
SIMC, in counts per mC, by the integrated beam cur- 
rent and then using the same filling procedure. 

As can be seen in Fig. [TSJ the resolution of the experi- 
ment, and therefore the width of the r\ peak, depends on 
cosO*. The integration limits, ji ow and j'high — the dashed 
lines in the figure, are also functions of cos9*. The depen- 
dence of the extracted cross-sections on these integration 
limits is accounted for in the next section. The exper- 
imental cross-section is then obtained, from the model 



cross-section at the centre of the bin a 



Mo- 



using 
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i data „i 

T data — Th ct MC- 



TV, 



(20) 



MC 



As with any measurement in which the events are 
histogrammed, each bin represents a mean quantity, 
weighted by the distribution of the events within that 
bin. In this experiment the cross-section changes rapidly 
and non-linear ly with W, especially going from thresh- 
old to maximum within just 50 MeV, and our W bins are 
rather large at 30 MeV. 

The bin centering in W was done implicitly during the 
cross-section extraction, under the assumption that the 
relativistic Breit-Wigner model and SIMC are accurate 
representations of the physics and detector response. If 
the simulation experimental acceptance model is accu- 
rate, then the kinematic distribution of simulated parti- 
cles in each bin will mimic the population of data events 
within that bin. So too if the physics model is good, 
then nonlinearities in the actual cross-section will be cor- 
rectly reproduced by the simulation. To the extent that 
both of these are true, the ratio of the data and Monte 
Carlo yields in each bin, A r data/-^MC> directly connects 
the number of detected particles with the input Monte 
Carlo model, uJj C . The bin centring is then done by 
evaluating the simulation input model at the bin centre. 

The bins in cos#* and </>* are quite small, and where 
there is full coverage, the extracted differential cross- 
sections are largely flat. It was decided not to attempt 
to incorporate nonlinear variation of the angular cross- 
section into the input model, and thus no implicit bin 
centering takes place. 

The results are not quoted at fixed Q 2 . Since the events 
in every bin have a Q 2 distribution, the cross-section re- 
sults are an average over the Q 2 distribution of the bin. 
The weighted average Q 2 of events in each bin (Q 2 ,^), is 
therefore quoted along with the extracted cross-section 
in Tables IXIII and IXIII1 In order to quote all the data at 
a single value of Q 2 , a model dependent correction would 
have to be applied to the data, which can be done at a 
later stage. 



E. Check of SOS acceptance 



incidence events were detected. At Ohms — 18° and 
Phms = 4.7 GeV/c, the acceptance for elastic ep events 
is best matched. At Ohms — 19.5° and Phms — 
4.5 GeV/c, the HMS in-plane angular acceptance reduces 
the SOS in-plane angular range to 44° to 47.5°. While 
for Ohms = 16.5° and P H ms = 4.5 GeV/c, the HMS in- 
plane angular acceptance reduces the SOS in-plane an- 
gular range to 49° to 51°. 

To extract measured elastic ep yields, the same data 
cuts listed in Table HVl were used with an additional cut of 
0.8 < W < 1.07 GeV to isolate elastic events. The data 
were also corrected for tracking efficiency, trigger inef- 
ficiency, computer and electronic deadtime. The same 
SIMC Monte Carlo was used with ep elastic cross section 
calculated using the electric and magnetic form factors 
from the fit of Bosted [H. At this Q 2 = 6.76 (GeV/c) 2 , 
the proton magnetic form factor is the dominant con- 
tribution to the elastic cross section and a conservative 
estimated error on the predicted cross section is 4%. 

8(%) 

20 18 16 14 12 10 8 6 4 2 

J I I I I I I U I L 



u 



IT 0.9 

>-■ 

"5 
| 0.8 

t3 
o 
£0.7 



~7 



mr 



t 



> HMS =18.0 andp HMS = 4.7GeV/c 
Ws = 19.5° and p„ MS = 4.5 GeV/c_ 
i 9„,„= 16.5° and p„= 4.5 GeV/c 



HMS 



'45 46 47 48 49 50 51 

6 C (deg) 

FIG. 19: (Colour online) Ratio of yield of elastic ep coin- 
cidence events to predicted yield from Monte Carlo (Yield 
Data/Yield MC) plotted versus 6 e for Osos - 47.5° and three 
different combinations of Ohms and Phms- The solid line is 
the average ratio = 0.95 ± 0.01, of all points between 6 e = 
45.5° to 49.5°. The corresponding value of electron 8 for a 
given 9 e is given by the upper x-axis. 



1. Coincident Elastic Scattering Cross-Section 

For the SOS central momentum and angle setting of 
SOS = 47.5° and P S os = 1-74 GeV/c, the scattered 
protons from elastic ep events will have a momentum 
of 4.44 GeV/c and angle of 18.3°. The elastic electrons 
cover an electron momentum range of 2.08 to 1.73 GeV/c 
and angular range of 44° to 51° which corresponds to 
a proton momentum range of 4.25 to 4.61 GeV/c and 
angular range of 19.8° to 17.0°. The Q 2 range is from 
6.4 to 7.1 (GeV/c) 2 . 

During the experiment, the HMS was set at three com- 
binations of Ohms and Phms at which elastic ep co- 



in Fig. Q33 the ratio of data yield to predicted Monte 
Carlo yield is plotted as a function of electron scattering 
angle for all three settings. Between scattered electron 
angle of 45.5° to 49.5°, the ratio is reasonably constant 
with an average value of 0.95 ± 0.01 which indicate good 
agreement with previous measurements. Below 45.5°, 
the agreement falls off sharply and above 49.5° the ratio 
jumps to an average of 1.08 which demonstrate problems 
in understanding the SOS acceptance in some areas. 

On the other hand, Figure Q13 shows that we are able 
to reproduce a well known quantity, the elastic scatter- 
ing cross section, to within a few percent using our two 
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spectrometer coincidence configuration, and thus we de- 
velop some confidence in the main result of the paper. 
The elastic events are in the SOS relative momentum 
range 10% < 5 < 20%, while the p(e, ep)r\ cross section 
is extracted in the range —20% < 6 < —5%, thus we can- 
not use this data to correct the rj cross section. For this 
reason a single arm comparison is best for checking the 
SOS acceptance. The ep coincidence comparison is use- 
ful as a check on the understanding of the experimental 
luminosity and efficiency corrections. 



2. Inclusive Elastic and Inelastic Cross-Section 

In order to determine how accurately the SIMC sim- 
ulation package models the acceptance of the SOS spec- 
trometer, we extracted single-arm elastic and inelastic 
cross-sections from hydrogen and compared them with 
a fit to previous data. This inclusive analysis had the 
same set of data runs, the same correction factors when- 
ever applicable, the same acceptance simulation code and 
the same electron identification cuts, as the coincidence 
analysis. 

In the inclusive case, corrections for the target end- 
caps were much larger than in the coincidence case, and 
an additional correction for pair-symmetric backgrounds 
was needed (up to 10% at the highest W). These were de- 
termined using interpolated positron cross sections mea- 
sured in a previous experiment [311 ] with the same target 
and beam energy, but slightly different scattering angles 
at 45, 55, and 70 degrees. This correction is negligible for 
the coincidence analysis due to the imposition of missing 
mass cuts. 

Another difference is that radiative corrections were 
done analytically, rather than in the Monte Carlo sim- 
ulation. For both elastic and inelastic scattering these 
were calculated using the formalism of Mo and Tsai [32| • 
For the required elastic scattering cross-section model, 
we used the form factor parametrisation of Bosted [3(| , 
while for the inelastic cross-section model we used the 
May 2007 fit of Christy and Bosted [H. 

To obtain final radiated cross-section for a proton tar- 
get, the cross-sections from the Al dummy target were 
subtracted with the appropriate scale factor to match 
the thickness of the endcaps. The small difference in ra- 
diative corrections between the endcaps and dummy was 
not taken into account. 

The W dependence of the extracted inelastic cross- 
section, taken from the central region of the SOS spec- 
trometer is plotted in Fig. [20] along with the Christy 
model. Generally, the ^-dependence is in quite good 
agreement with Christy fit, especially for 1.5 < W < 1.6 
GeV, which is the main focus of this paper. Based on this 
analysis, a systematic uncertainty of 3% was assigned to 
the acceptance of the SOS spectrometer. 
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FIG. 20: (Colour online) Inclusive inelastic differential cross- 
sections as measured by the SOS spectrometer centered at 
47.5 and 70 degrees, as a function of W , with the angular cut 
— 30 < dy/dz < 30 mr. The curves arc from a fit to world 
data El. 



F. Systematic Error Analysis 

Depending on the source of error, one of two differ- 
ent methods was used to account for it. Those errors 
that were independent of the kinematic variables of the 
extracted data, W, cos8* and (f>* were treated globally 
and applied to the data overall. The sources of this kind 
of error are summarised in Tabic I VIII 

The pion contamination through the PID cut for elec- 
trons was calculated by Villano 20], using the data of 
this experiment, to be 1.6%. Most of these pions are 
from random coincidences and are effectively removed 
by the coincidence time cut — an analysis for the F T ex- 
periment pjj [21j shows the residual contamination to be 
about 0.1%. The systemstic error in the target density 
and charge measurement were also determined by the 
F w analysis [lj| |2l| . The error in the HMS acceptance 
is the quadrature sum of the 0.5% point-to-point error 
and 0.8% normalisation error determined by Christy [331 ] . 
The overall error of 4.2%, calculated as a sum in quadra- 
ture, is dominated by the uncertainty in the SOS accep- 
tance. 

If a source of error was expected to be dependent on 
kinematics, then it was treated on a bin-by-bin basis. 
The Monte Carlo simulation was run with altered param- 
eters to mimic the uncertainty, and the subsequent anal- 
ysis was done to compare to the extracted cross-section 
and quantify the effect bin-by-bin. The parameters that 
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Uncertainty Reference 



Parameter 

SOS acceptance 3.0% " Sec. IIV E 21 

Radiative Corrections 2.0% Sec. IIV C 51 

Trigger efficiency 1.4% Sec. IIV A II 

Proton absorption 1.0% Sec. IIV All 

HMS acceptance 1.0% Ref. [33] 

Target density 0.6% Ref. [19] 

Charge measurement 0.5% Ref. [191 

Electron PIP cut 0.1% Ref. [19, 20] 

TOTAL (quadrature sum) 4.2% 



TABLE VII: The sources of global systematic error and their 
estimated sizes. 



were altered, listed in Table IVlTTI were those considered 
imprecisely known or that affect the determination of the 
cross-section. The best choice set of parameters were 
used for the standard analysis from which the final differ- 
ential cross-section was calculated. Each parameter was 
then varied and the complete analysis repeated, up to 
the point of attaining the differential cross-section. The 
parameters were not varied together, as would be done in 
a fit, since it was assumed that to first order they acted 
independently and thus the prohibitive extra effort was 
unnecessary. 

The drift chamber resolutions, r^c, for the HMS and 
SOS spectrometers were calibrated as described in Sec- 
tion IIV A 21 In order to completely account for any er- 
ror, these parameters were arbitrarily increased by 10% 
for the variation procedure. The exact position of the 
target in the beam direction, z tar g, was only known to 
within 3 mm. For the standard analysis, the middle po- 
sition of this uncertainty window, an offset of 1.5 mm 
from the nominal centre, was chosen. The variation used 
for this parameter was the maximum possible extent of 
the motion, 1.5 mm in cither direction. 

The SOS spectrometer was found to be somewhat out- 
of-plane, but the exact amount is uncertain. A survey 
of the hall produced a value of x' sos = 2.62 mr, which 
was used in this extraction, while an analysis of ep co- 
incidence data by the F„ experiment [ijj, [2l| yielded 
x' sos = 3.2 mr. The spectrometer offset was thus var- 
ied in both directions, to 1.5 mr and 3.5 mr, for the 
systematic analysis. 

The cut on missing mass squared ra 2 ., is described in 
Sec. IIV B 31 The effect of this cut was taken into account 
by including it as one of the parameters varied in the 
systematic analysis. The variation chosen was to widen 
this cut on both ends by 0.1 GeV 2 and then subsequently 
to narrow it by the same amount. 

If Xi was the value of the differential cross-section in bin 
i for the standard analysis and y\ was for the analysis of 
a certain variation v, then the systematic error for that 
variation in that bin was taken as half the difference, 
8? = \xi-yX\/2. 

For the purposes of conveying the size of each of the 
systematic errors in Table IVIIII a measure of the average 



size (S v ) is used. This is the mean systematic error for all 
bins, weighted by the statistical error of the measurement 
in each bin 



<<n 






where Ui is the statistical error of the differential cross- 
section in bin i. 



Parameter 




Pstd 




pvar 


in 


HMS mc 


(mm) 


0.57 




0.66 


3.1% 


SOS r DC 


(mm) 


0.35 




0.39 


3.7% 


z S os offset 


(mr) 


2.62 




1.5 


3.1% 






2.62 




3.5 


2.8% 


Ztarg Offset 


(mm) 


1.5 




0.0 


3.0% 






1.5 




3.0 


2.9% 


m 2 , cut 


(GeV 2 ) /(cos6»* 


')/, 


max+0.1 
Tiin-0.1 


3.5% 






f(cos6' 


") /, 


max — 0.1 
nin+0.1 


2.8% 



TABLE VIII: The various sources of kinematic dependent sys- 
tematic errors considered in the analysis, the standard sim- 
ulation values Pstd, the systematic variation p var , and the 
weighted mean systematic error for all bins, (S v ). 

The total bin i systematic error, <5* ot , was determined 
by adding in quadrature the systematic error for each 
variation, 8%, and the global systematic errors, S g i , to 

give <5* ot = 



E^) 2 + Mo- 



V. RESULTS 
A. Differential Cross-Section p(e, e'p)rj 

The differential cross-sections for the centre-of-mass 
scattering angles of the r\ are extracted in the bins de- 
scribed in Section IIV A 51 with large W bins to allow 
more angular bins. Figure HH shows these data for the 
lower-Q 2 setting. The diminishing experimental accep- 
tance as W increases, especially in out-of-plane <fi* bins, 
is evident. As seen in previous data @, [f| [1> 03 > a domi- 
nant isotropic, or 5-wave, component is seen at W from 
threshold to the Sn resonance peak. 

Equation ([8]) is the parametrisation of the virtual pho- 
ton cross-section in terms of its angular dependence. 
The extracted differential cross-section was fitted with 
Eq. © , for the lower W bins where there is sufficient an- 
gular acceptance for a fit, and is plotted in Fig. [21] The 
parameters extracted from the fit are plotted in Fig. |2"21 
and listed in Table ITXl Using the results of the fit, the 
anisotropy in the threshold to resonance region is shown 
to be at most about 15% for the lower-Q 2 setting. 

The results of this fit can be compared to similar stud- 
ies of the angular dependence of rj production data. The 
recent CLAS data [6] was also fit with Eq. ([5]). The term 
linear in cos9* shows definite structure at all measured 
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FIG. 21: (Colour online) Extracted ep — > epr; differential cross-sections for the lower-Q setting. The solid (blue) curve is a 
fit of Eq. (JSj to each W bin. The dashed curve is the eta-maid [3J isobar model for ry-electroproduction from the nucleon 
at Q 2 = 5 GeV 2 , projected to the appropriate Q 2 for each W bin by the factor (5 GeV 2 /Q 2 (W / )) 3 . The inner error bars are 
statistical and the outer error bars, the quadrature sum of the statistical and systematic errors. 
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W = 1500 MeV W = 1530 MeV W = 1560 MeV W = 1590 MgV W = 1625 MeV 


A 


63.34 ± 1.59 


71.25 ± 1.53 


47.75 ± 1.41 


30.65 ± 1.50 


19.84 ± 3.42 


B 


2.70 ± 1.89 


10.95 ± 1.88 


7.54 ± 1.79 


-1.74 ± 2.18 


-5.53 ± 7.81 


C 


7.27 ± 3.71 


8.17 ± 3.62 


21.37 ± 3.37 


14.63 ± 3.66 


7.88 ± 6.58 


D 


-2.53 ± 1.89 


-0.57 ± 1.83 


0.10 ± 1.66 


1.67 ± 1.83 


2.80 ± 5.06 


E 


-3.83 ± 4.16 


-3.67 ± 4.04 


4.11 ± 3.59 


3.47 ± 3.83 


4.01 ± 8.16 


F 


4.86 ± 2.04 


8.17 ± 1.96 


5.06 ± 1.78 


5.73 ± 1.86 


4.53 ± 2.65 



TABLE IX: The extracted angular parameters from a fit of Eq. (jHJ to the lower-Q 2 extracted differential cross-section. 
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FIG. 22: (Colour online) Extracted parameters from fits of 
Eq. ((8} to the lower-Q 2 differential cross-section, shown as 
curves in Fig. [21] 



FIG. 23: (Colour online) The result of fits to the differen- 
tial cross-section, plotted as the ratio of the linear cosf?^ term 
to the isotropic component, for the present work and other 
?7-electroproduction data |2, la Hg|> The black dotted line is 
drawn at W — 1.535 GeV, the nominal mass for the Sn res- 



Q 2 . It was observed that as W increases above where 
the 5n(1535) is expected to be dominant, the cos6>* de- 
pendence changes dramatically. At W — 1.66 GeV it 
decreases monotonically with cos#* but by W — 1.72 
GeV the forward backward-asymmetry is reversed. Pre- 
vious experiments, at photoproduction [35| and at higher 
Q 2 @, have shown the same structure in the W depen- 
dence of B, with B J A appearing to be roughly indepen- 
dent of Q 2 up to Q 2 = 2.5 GeV 2 @. 

The quantity B/A for the present work and previously 
published data 0, la, Q is plotted in Fig. [23] Due to 
diminishing angular acceptance the present work does 
not extend above W ~ 1.65 GeV where the ratio reaches 
its minimum and begins to make a rapid change from 
negative to positive. For W near the S\\ resonance mass 
(black dotted line in figure), the B/A structure shows 
some difference between the CLAS data [fj| which remains 
negative and data from the present work and others [2, 
[5J which do go positive, but the trend is the same and 
continues to be approximately independent of Q 2 up to 
-5.8 GeV 2 . 

The higher-Q 2 setting data were not amenable to the 
full angular fit, as can be seen in Fig. [Ml so the fit func- 
tion was truncated to da/dVl* = A + B cos8* and fitted 



to the data. There is large uncertainty on the extraction 
of B/A for these data, and the results are consistent with 
no structure, as can be seen in Fig. [23] 

Denizli et al. 6] show that the rapid change in sign 
of B could be due to a P wave resonance at W ~ 1.7 
GeV. Specifically, a simple resonance model incorporat- 
ing the Pn(1710) could describe their data, but they do 
acknowledge that the Pi3(1720) is also a candidate. The 
approximate Q 2 independence of the magnitude of this 
feature would imply that such a P wave falls similarly 
slowly with Q 2 as the Sii(1535). 

As can be seen in Fig. [211 the cos 2 #* term in the angu- 
lar fit to the lower-Q 2 data is also quite significant for W 
above the resonance mass. In this case, the agreement 
with [y] is not good, as can be seen in Fig. [25j This dis- 
agreement can also clearly be seen qualitatively in Fig.[2Tl 
where the eta-maid [3J| curves are concave down while 
the new data are concave up. 
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FIG. 24: (Colour online) Extracted ep — > ep?7 differential 
cross-sections for the higher-Q 2 setting. The (blue) solid 
curve is a fit to the data of the form da/dQ* — A + B cos9* . 
The dashed curve is the eta-maid model [H| at Q 2 = 5 GeV 2 , 
projected to the appropriate Q 2 for each W bin by the factor 
(5 GeV 2 / Q 2 (W)) 3 . The inner error bars are statistical and 
the outer error bars, the quadrature sum of the statistical and 
systematic errors. 



B. Total Cross-Section p(e, e!p)r\ 




1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 
W [GeV] 



FIG. 25: (Colour online) The ratio of the quadratic 
cos 2 6* term to the isotropic component for fits to the r\- 
electroproduction differential cross-section, for the present 
work and other data @, [f|- Symbols are the same as Fig. I23I 



(Q 2 ) [*¥■] 


W [GeV 


a [nb] 


Y.%1^ 


5.802 


1.50 


831.9 ± 19.7 


100.0% 


5.764 


1.53 


926.5 ± 20.9 


100.0% 


5.704 
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100.0% 


5.636 


1.59 


461.0 ± 15.0 


93.8% 


5.554 


1.62 
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45.8% 


5.456 


1.67 


247.8 ± 14.7 


29.2% 


5.353 


1.71 


239.0 ± 17.1 


18.8% 


5.248 


1.75 


175.6 ± 17.8 


16.7% 


5.136 


1.78 


160.3 ± 31.7 
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5.022 


1.83 


162.9 ± 27.2 


7.3% 


7.064 
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482.1 ± 33.3 


43.3% 


7.011 


1.53 
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36.7% 
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33.3% 


6.857 


1.59 


282.6 ± 25.0 


23.3% 


6.746 


1.64 


228.9 ± 26.9 


16.7% 


6.602 


1.69 


168.1 ± 37.5 


10.0% 


6.462 


1.74 


230.3 ± 60.3 


3.3% 



The total cross-section was determined from the differ- 
ential cross-section in two ways. Firstly, the total cross- 
section was obtained by taking the weighted mean of the 
differential cross-section in each W bin and multiplying it 
by 47r, where the uncertainty in the mean is the quadra- 
ture sum of the statistical and systematic errors from all 
the bins. In W bins where there is full coverage, this 
is equivalent to integrating the differential cross-section. 
The total cross-section found using this method is listed 
in Table IXl along with the weighted average Q 2 in each 
W bin and the percentage of the 47r cm. angular range 
accepted in each W bin. Secondly, the fitted angular de- 
pendence, Eq. © with parameters given in Table IJ"Xl was 
integrated in each W bin. Here, the uncertainty was de- 
termined by fixing each of the six parameters to the high 
and low one-sigma Minuit fit values and then fitting the 
remaining five parameters and determining the integral. 



TABLE X: Table of the total cross section, determined from 
the weighted average of extracted differential cross section. 
The weighted average Q 2 and the percentage of angular cov- 
erage for each W bin are also indicated. The errors are sta- 
tistical and systematic added in quadrature, and do not take 
into account the angular acceptance. 



The maximum and minimum values of the integral so de- 
termined were used to estimate the error. This second 
procedure couldn't be applied to the higher-Q 2 setting 
because the sparsity of the data precluded the fitting of 
the full angular dependence. The total cross sections de- 
termined in this way for each of the settings still have a 
Q 2 which varies with W. 

The value of <jr was obtained by fitting a relativis- 
ts Breit-Wigner to the total cross-section and evaluat- 
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ing it at the resonance mass. The Breit-Wigner is given 
by Eq. [T7| and described in Section IIVC II and the non- 
resonant background is modeled as A ni y/W — Wthr + 
B n i(W — Wthr)- During the fit, the mean Q 2 value for 
that W bin was used. Due to strong correlations between 
the parameters, especially b^ and Wr, the branching frac- 
tion to 77 was fixed at b^ = 0.5 for the fits. The uncer- 
tainty in ur was estimated by individually fixing each of 
the Breit-Wigner parameters Wr and Tr to their Minuit 
uncertainties, redoing the fit and reevaluating ur. The 
maximum and minimum values of or so determined were 
used to estimate the error. 

This method worked well for the lowerQ 2 data, with 
good agreement of a single Breit-Wigner to the data. For 
the averaged differential cross-section a small background 
contribution, less than 0.5%, was admitted under the 
resonance peak, while the fit to the integrated angular 
dependence model didn't admit any background contri- 
bution. The higher-Q 2 data were amenable to such a 
fit since the large error bars and poor angular coverage 
make the parameters unreliable. For this reason, a si- 
multaneous fit to both settings was thus done, yielding 
a single set of resonance parameters. The background 
was constrained to have the same Q 2 dependence as the 
data, essentially requiring it to have the same relative 
size. Figure [26] shows the results of this fit, which are 
listed in Table IXj along with the results of the fits to the 
lower-Q 2 data. The shape of the fitted function is domi- 
nated by the lower-Q 2 data, with a background of 1.2% 
at the resonance mass. The values from this simultane- 
ous fit are used in the further analysis. 
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FIG. 26: (Colour online) A simultaneous fit to the lower-Q 2 
and higher-Q 2 data of the sum (solid line) of a relativistic 
Breit-Wigner (long dash) and non-resonant background term 
(short dashed line). The data are the total cross section de- 
termined from 4n{d<j/dQ*). The background was constrained 
as described in the text. 



Both the simultaneous and the individual fits were re- 
peated for b v = 0.45 and 0.55. The results of these addi- 



tional fits are plotted as correlation contours in Fig. [27] 
The or extracted from each of these additional fits was 
at all times well within the error quoted in Table IXfl ft 



can be seen that there are correlations between b n and 
Wr and also between Wr and Tr. The resonance pa- 
rameters from the simultaneous fit are dominated by the 
lower-Q 2 data, as expected. 
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FIG. 27: (Colour online) Plot of the 1-sigma contours from 
the various Breit-Wigner fits to the data. 



Helicity Amplitude A\/2 for the 511(1535) 
Resonance 



The amplitude A 1 / 2 is determined from the total 
cross-section at the 5*11(1535) resonance mass ur, by 
Eq. (fl~2"]) , which assumes A 1 / 2 3> S±/ 2 - Using the <tr 
values obtained from the Breit-Wigner fit to the to- 
tal cross-section, and those obtained in previous exper- 
iments @, [fl, Q, A-^/2 is determined consistently for all 
data with T R = 150 MeV, b n = 0.55 and Wr = 1535 
MeV, chosen to coincide with those used previously @,@|- 
The uncertainties in A\ / 2 do not include uncertainties in 
Wr, b n or Tr. 

Table IXII summarises the parameters from the Breit- 
Wigner fit, the extracted total cross-section at the res- 
onance mass, ur, and the extracted helicity ampli- 
tude, Ax/2- As can be seen in Fig. [25] the values of 
A1/2 determined in this work significantly extend the 
Q 2 ran ge o f the world's data. The curves in the fig- 
ure [3a, [33, [33, [33, EO] show a huge variation in the pre- 
dicted values of Ax/2- 

The magnetic form-factor of the proton G P M as pub- 
lished by Arnold et al. [4JJ demonstrates clear scaling 
behaviour. Naive dimension counting in pQCD predicts 
a falloff of 1/Q 4 and the quantity of Q A G V M reaches a 
broad maximum at about Q 2 ~ 8 GeV 2 and then de- 
creases in a gentle logarithm due to the running of the 
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Q\Wr) 
[GeV 2 /c 2 ] 



J'dcr/df7*| modc i indiv. 
47r(da/df2*) 



W R 

[GeV] 



r R 

[GeV] 



OR 
\pb\ 



Al/2 

fxlO" 3 GeV 5/2 l 



5.79 1.523 ± 0.001 0.125 ± 0.003 0.976 ± 0.007 23.62 ± 0.09 
5.79 1.523 ±0.001 0.128 ±0.010 0.977 ±0.024 23.63 ± 0.30 



Fig. 1261 5/79 1.522 ±0.001 0.128 ±0.009 0.943 ± 0.015 23.22 ±0.18 
Fig. HD 7.04 0.553 ±0.020 17.79 ± 0.33 



47r(dcr/dfi*) 
A-K(do/dQ,*) 



TABLE XI: Parameters extracted from a relativistic Breit-Wigner fit to the data. The values for A\/2{Q 2 ) are determined 
from an, assuming A\/ 2 2> Sj/2, with parameters Wr — 1.53 GeV, Tr = 150 MeV and b v — 0.55. 
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FIG. 28: (Colour online) Values for A V| 
from or for the present and other data [3 ; 
with W R = 1.53 GeV J" fl = 150 MeV and b 
curves are from Refs. Ho, 013. 138, 39, 401. 



2 ) determined 

(consistently 

= 0.55). The 



FIG. 29: (Colour online) The Q 2 dependence of Q 3 A 1/2 for 
^-production. Scaling in this quantity appears to begin at a 
photon momentum transfer of Q 2 ~ 5 GeV 2 . The dashed lines 
are a high Q 2 , pQCD calculation from Carlson and Poor [J] 
using three different nucleon distribution amplitudes. 



strong coupling constant a s . The same arguments pre- 
dict that the helicity amplitude for the 5n(1535) de- 
creases with 1/Q 3 . Figurel2Ulis a plot of Q 3 Ai/ 2l showing 
that the quantity Q 3 Ai/ 2 appears to begin flattening at 
a photon momentum transfer broadly within the range 
of this work, Q 2 ~ 5 — 7 GeV 2 , a possible signal of the 
onset of pQCD scaling. A pQCD calculation by Carl- 
son and Poor [4(], of the magnitude of this quantity, is 
plotted and is a factor of ~ 3 smaller than the data. It 
has also been pointed out that such scaling may have a 
non-perturbative explanation [42|, |43j . 



In order to compare the behaviour of A 1 / 2 with the 
approach of G P M to scaling, the quantity Q 3 A 1 / 2 /Q 4 G p M 
is plotted in Fig. [301 The form of G P M is taken from the fit 
by Bosted [3fj]. The figure shows that the two quantities 
don't have the same form at low Q 2 , and the data doesn't 
go high enough in Q 2 to know whether the two quantities 
begin behaving equivalently. 
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FIG. 30: (Colour online) The Q 2 dependence of 
Q 3 A 1/2 {Su)/Q 4 G p M for ^-production. 
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APPENDIX: TABLES OF DIFFERENTIAL 
CROSS-SECTIONS 



VI. CONCLUSIONS 

We have presented the results of a precise, high statis- 
tics measurement of the differential cross-section for the 
ep — > e'prj exclusive process. This is done at the highest 
momentum transfer to date, namely, Q 2 = 5.8 and 7.0 
(GeV/c) 2 at the Sn resonance mass, which is a signif- 
icant extension from the previous highest at Q 2 = 3.6 
GeV 2 . Data were obtained from threshold to W = 1.8 
GeV, the 5n(1535) dominating the channel as expected. 
In the region from threshold to the S l n(1535) resonance 
mass, the differential cross-section is largely isotropic — 
consistent with previous measurements. 

The interference phenomenon in the linear cos6>; term TABLE XII: Lower-Q 2 extracted differential cross-section, 

at W of the 5n(1535) resonance mass, seen in lower Q 2 
and photoproduction data is observed here with simi- 
lar strength. The present data doesn't have sufficient 
angular coverage at W ~ 1.7 GeV to comment mean- 
ingfully on the strong presence of a P wave resonance 
there. The curvature in the cos9* dependence of the dif- 
ferential cross-section is opposite to that of the data at 
lower-Q 2 - The helicity-conserving transition amplitude 
Ai/2, is extracted from the data assuming no longitu- 
dinal component (A l / 2 3> <Si/ 2 ). The Q 2 dependence of 
Q 3 Ai/2 seems to be flattening, consistent with the pQCD 
prediction, although the range of Q 2 is too small to verify 
the exact dependence. Even if the data scale as predicted 
by pQCD, that is not conclusive evidence for the onset 
of pQCD. 

On the theoretical front, the differential cross-section 
will be incorporated into mult i- channel, multi-resonance 
models, such as those by the MAID and EBAC groups, 
which should maximize the physics impact coming from 
these data. Also, the inability for any one calculation 
to adequately describe the Q 2 dependence of A 1 / 2 leaves 
much to be done in understanding the structure of the 
S'ii(1535). On the experimental front, more data are 
required to further address the questions in this paper. 

It would be nice to fill the data gap in the region be- 
tween Q 2 ~ 4 and 5.8 GeV 2 to analyse the apparent 
change of differential cross-section shape. Extending the 
data to Q 2 much higher than 7 GeV 2 will complete the 
study of the transition to hard-scale scattering. Obtain- 
ing LT separated data at high Q 2 will enable checking 
of the assumption, made in this work and in the litera- 
ture, that the longitudinal component is negligible. The 
planned upgrade of the Jefferson Lab accelerator, to en- 
ergies as high as 11 GeV, will allow exclusive r/ electro- 
production data to be obtained to Q 2 ~ 14 GeV 2 , and 
LT separations at least to the Q 2 of this experiment. 
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1.530 0.250 247.5 5.77 0.423 1.27 70.3 9.1 3.3 

1.530 0.250 292.5 5.77 0.423 1.31 67.2 9.7 3.4 

1.530 0.250 337.5 5.76 0.425 1.40 81.6 10.4 3.9 

1.530 0.417 22.5 5.76 0.425 1.37 72.1 9.4 4.5 
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1.530 0.417 247.5 5.77 0.422 1.29 75.0 8.5 4.4 



27 



W cosO* 

[GcV] 



[deg. 



(QlJ <^> 

rGeV 2 i 



rad . 
corr. 



d*cr 
dS2* 

[ — ] 

L sr -I 



Cstat Osyst 
r nb l r nb "i 



W cosi9* 

[GeV] 



[deg 



1 <Qbin) <£> 
GcV 2 i 



rad . 
corr. 



1.530 0.417 292.5 5.77 0.423 1.32 84.1 9.7 5.0 

1.530 0.417 337.5 5.77 0.424 1.38 94.7 10.3 4.8 

1.530 0.583 22.5 5.77 0.423 1.36 77.4 11.2 5.6 

1.530 0.583 67.5 5.77 0.423 1.27 81.7 10.4 4.7 

1.530 0.583 112.5 5.77 0.422 1.28 63.0 8.8 4.7 

1.530 0.583 157.5 5.77 0.422 1.30 84.8 10.5 4.6 

1.530 0.583 202.5 5.77 0.423 1.30 95.2 11.1 4.6 

1.530 0.583 247.5 5.77 0.423 1.29 94.1 10.2 5.4 

1.530 0.583 292.5 5.76 0.424 1.32 84.2 10.8 6.0 

1.530 0.583 337.5 5.77 0.423 1.40 72.9 10.5 3.7 

1.530 0.750 22.5 5.79 0.421 1.33 87.7 12.3 4.7 

1.530 0.750 67.5 5.77 0.423 1.31 73.2 11.0 4.4 

1.530 0.750 112.5 5.77 0.423 1.28 75.5 10.6 3.7 

1.530 0.750 157.5 5.77 0.422 1.29 93.1 11.7 7.2 

1.530 0.750 202.5 5.77 0.423 1.32 77.5 10.8 4.3 

1.530 0.750 247.5 5.77 0.423 1.32 82.0 10.1 4.6 

1.530 0.750 292.5 5.77 0.423 1.32 68.4 10.8 5.0 

1.530 0.750 337.5 5.79 0.421 1.35 78.4 11.3 6.3 

1.530 0.917 22.5 5.78 0.423 1.31 94.7 12.6 7.5 

1.530 0.917 67.5 5.76 0.425 1.31 86.5 11.8 5.7 

1.530 0.917 112.5 5.77 0.424 1.31 104.3 11.6 6.2 

1.530 0.917 157.5 5.77 0.423 1.30 77.4 10.5 5.2 

1.530 0.917 202.5 5.77 0.423 1.31 106.5 11.6 5.9 

1.530 0.917 247.5 5.77 0.423 1.31 102.1 11.5 5.1 

1.530 0.917 292.5 5.77 0.424 1.31 77.8 11.2 4.6 

1.530 0.917 337.5 5.77 0.423 1.31 106.0 12.8 7.4 

1.560 -0.917 22.5 5.70 0.421 1.38 67.3 7.1 4.7 

1.560 -0.917 67.5 5.70 0.421 1.36 62.2 7.0 3.7 

1.560 -0.917 112.5 5.70 0.420 1.35 62.5 7.0 4.4 

1.560 -0.917 157.5 5.69 0.422 1.33 67.8 7.0 4.4 

1.560 -0.917 202.5 5.69 0.422 1.32 55.5 6.2 3.2 

1.560 -0.917 247.5 5.70 0.421 1.34 59.5 6.2 4.3 

1.560 -0.917 292.5 5.69 0.422 1.36 54.1 6.1 3.5 

1.560 -0.917 337.5 5.70 0.420 1.38 55.7 6.6 3.2 

1.560 -0.750 22.5 5.70 0.422 1.34 60.3 6.6 2.9 

1.560 -0.750 67.5 5.69 0.422 1.32 53.6 6.3 3.8 

1.560 -0.750 112.5 5.69 0.423 1.28 50.0 5.9 4.6 

1.560 -0.750 157.5 5.68 0.423 1.25 54.2 6.1 3.8 

1.560 -0.750 202.5 5.69 0.422 1.24 50.9 5.6 3.1 

1.560 -0.750 247.5 5.69 0.422 1.31 54.6 5.9 2.8 

1.560 -0.750 292.5 5.70 0.420 1.32 50.1 5.8 3.1 

1.560 -0.750 337.5 5.70 0.421 1.35 52.9 6.1 2.8 

1.560 -0.583 22.5 5.69 0.423 1.39 49.3 6.5 2.6 

1.560 -0.583 67.5 5.71 0.420 1.34 46.8 6.9 2.9 

1.560 -0.583 112.5 5.70 0.422 1.24 44.9 6.6 3.9 

1.560 -0.583 157.5 5.70 0.421 1.27 49.4 6.1 3.4 

1.560 -0.583 202.5 5.69 0.421 1.24 54.0 6.0 2.8 

1.560 -0.583 247.5 5.70 0.421 1.26 42.0 5.8 2.5 

1.560 -0.583 292.5 5.71 0.419 1.30 55.6 6.8 2.8 

1.560 -0.583 337.5 5.69 0.422 1.36 52.5 6.4 2.8 

1.560 -0.417 22.5 5.68 0.425 1.33 45.1 7.4 2.8 

1.560 -0.417 67.5 5.72 0.419 1.26 49.9 7.7 2.9 

1.560 -0.417 112.5 5.70 0.421 1.23 43.7 7.0 2.7 

1.560 -0.417 157.5 5.69 0.421 1.26 53.6 6.9 3.0 

1.560 -0.417 202.5 5.70 0.420 1.24 57.5 6.8 3.7 

1.560 -0.417 247.5 5.69 0.422 1.17 38.6 6.3 3.0 

1.560 -0.417 292.5 5.72 0.419 1.24 46.3 6.9 3.1 

1.560 -0.417 337.5 5.68 0.425 1.30 52.6 7.5 2.9 

1.560 -0.250 22.5 5.68 0.424 1.28 54.9 9.1 2.9 

1.560 -0.250 67.5 5.71 0.419 1.20 42.9 7.5 3.1 

1.560 -0.250 112.5 5.70 0.420 1.16 40.5 7.2 2.1 



d*cr 
dS2* 

[ — ] 

I sr -I 



y stat ^syst 
r nb 1 r nb i 



1.560 -0.250 157.5 5.69 0.421 1.22 56.5 7.6 2.7 

1.560 -0.250 202.5 5.70 0.420 1.23 59.6 7.1 4.3 

1.560 -0.250 247.5 5.71 0.419 1.16 51.1 7.2 2.9 

1.560 -0.250 292.5 5.72 0.419 1.20 40.8 6.9 4.5 

1.560 -0.250 337.5 5.68 0.425 1.30 41.6 8.1 3.1 

1.560 -0.083 22.5 5.70 0.422 1.27 54.3 10.3 3.6 

1.560 -0.083 67.5 5.73 0.417 1.15 32.9 8.1 2.2 

1.560 -0.083 112.5 5.72 0.419 1.06 38.8 7.5 4.1 

1.560 -0.083 157.5 5.71 0.419 1.21 64.5 8.7 4.1 

1.560 -0.083 202.5 5.70 0.420 1.19 57.9 7.9 3.3 

1.560 -0.083 247.5 5.72 0.420 1.07 70.5 9.5 3.9 

1.560 -0.083 292.5 5.72 0.420 1.17 51.9 8.8 3.4 

1.560 -0.083 337.5 5.69 0.423 1.26 41.8 8.8 2.9 

1.560 0.083 22.5 5.68 0.426 1.24 61.5 10.1 3.2 

1.560 0.083 67.5 5.72 0.420 1.10 48.4 10.5 2.8 

1.560 0.083 112.5 5.72 0.419 1.00 50.7 10.3 5.6 

1.560 0.083 157.5 5.69 0.421 1.19 56.0 9.2 3.3 

1.560 0.083 202.5 5.70 0.420 1.19 60.6 8.7 3.3 

1.560 0.083 247.5 5.74 0.418 1.02 44.3 8.9 3.3 

1.560 0.083 292.5 5.73 0.419 1.11 54.5 10.2 3.2 

1.560 0.083 337.5 5.69 0.424 1.27 56.8 9.7 3.4 

1.560 0.250 22.5 5.69 0.424 1.23 44.5 9.0 2.5 

1.560 0.250 67.5 5.73 0.419 1.06 30.8 9.7 2.0 

1.560 0.250 112.5 5.74 0.417 0.98 47.0 9.9 2.7 

1.560 0.250 157.5 5.70 0.420 1.14 43.3 8.8 2.6 

1.560 0.250 202.5 5.70 0.421 1.16 48.1 8.4 3.2 

1.560 0.250 247.5 5.74 0.417 0.99 51.9 9.4 2.7 

1.560 0.250 292.5 5.74 0.418 1.07 55.8 10.8 2.8 

1.560 0.250 337.5 5.69 0.424 1.25 47.9 9.1 2.2 

1.560 0.417 22.5 5.69 0.424 1.20 49.5 9.5 2.6 

1.560 0.417 67.5 5.73 0.418 1.05 49.8 10.7 3.1 

1.560 0.417 112.5 5.73 0.418 1.00 48.4 9.4 4.8 

1.560 0.417 157.5 5.70 0.420 1.13 51.8 9.3 3.1 

1.560 0.417 202.5 5.70 0.420 1.13 50.0 8.8 4.8 

1.560 0.417 247.5 5.73 0.418 1.01 53.0 9.4 4.1 

1.560 0.417 292.5 5.73 0.418 1.10 55.3 10.3 3.6 

1.560 0.417 337.5 5.69 0.424 1.25 64.1 9.7 3.7 

1.560 0.583 22.5 5.71 0.420 1.17 65.6 11.9 4.3 

1.560 0.583 67.5 5.74 0.416 1.05 61.8 11.8 4.0 

1.560 0.583 112.5 5.72 0.418 1.03 53.2 10.8 2.9 

1.560 0.583 157.5 5.71 0.419 1.11 56.0 11.3 4.5 

1.560 0.583 202.5 5.71 0.419 1.13 53.0 10.3 2.5 

1.560 0.583 247.5 5.72 0.418 1.05 44.3 10.1 3.8 

1.560 0.583 292.5 5.74 0.417 1.08 54.8 11.3 3.4 

1.560 0.583 337.5 5.71 0.421 1.21 63.4 11.0 3.3 

1.560 0.750 22.5 5.71 0.421 1.12 75.8 14.4 5.6 

1.560 0.750 67.5 5.70 0.422 1.09 67.5 13.0 4.2 

1.560 0.750 112.5 5.70 0.421 1.08 62.0 11.2 4.2 

1.560 0.750 157.5 5.71 0.420 1.14 60.4 10.6 2.8 

1.560 0.750 202.5 5.71 0.420 1.15 67.9 10.3 3.5 

1.560 0.750 247.5 5.72 0.419 1.12 67.2 10.1 5.5 

1.560 0.750 292.5 5.71 0.420 1.13 65.2 11.6 4.2 

1.560 0.750 337.5 5.73 0.419 1.12 86.2 13.5 4.0 

1.560 0.917 22.5 5.69 0.424 1.19 88.8 14.5 4.3 

1.560 0.917 67.5 5.69 0.424 1.17 87.7 13.1 4.9 

1.560 0.917 112.5 5.70 0.422 1.17 65.3 11.6 4.4 

1.560 0.917 157.5 5.70 0.422 1.18 72.7 11.8 3.6 

1.560 0.917 202.5 5.70 0.421 1.17 67.3 11.2 3.8 

1.560 0.917 247.5 5.70 0.422 1.19 64.9 11.1 4.4 

1.560 0.917 292.5 5.69 0.424 1.18 81.9 12.9 4.5 

1.560 0.917 337.5 5.69 0.425 1.21 77.9 14.3 4.5 
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W cosO* 

[GeV] 



[deg. 



(QlJ <^> 

rGeV 2 i 



rad . 
corr. 



d*cr 
dS2* 

[ — ] 

L sr -I 



fstat Osyst 
r nb 1 r nb i 



W cosi9* 

[GeV] 



[deg 



V (Qbin) <£> 
GcV 2 i 



rad . 
corr. 



1.590 -0.917 22.5 5.64 0.417 1.41 61.2 6.6 3.9 

1.590 -0.917 67.5 5.63 0.418 1.37 46.9 6.2 2.6 

1.590 -0.917 112.5 5.63 0.418 1.32 63.8 7.2 5.0 

1.590 -0.917 157.5 5.63 0.419 1.33 37.4 5.8 2.9 

1.590 -0.917 202.5 5.63 0.419 1.32 45.1 5.8 3.0 

1.590 -0.917 247.5 5.63 0.418 1.29 41.3 5.5 2.3 

1.590 -0.917 292.5 5.63 0.420 1.34 42.4 5.6 2.7 

1.590 -0.917 337.5 5.64 0.417 1.38 42.5 5.5 2.6 

1.590 -0.750 22.5 5.62 0.421 1.31 49.0 6.2 2.8 

1.590 -0.750 67.5 5.64 0.418 1.29 47.8 6.3 3.2 

1.590 -0.750 112.5 5.63 0.418 1.24 49.0 6.3 2.9 

1.590 -0.750 157.5 5.63 0.419 1.21 47.4 5.8 2.6 

1.590 -0.750 202.5 5.63 0.419 1.22 40.4 5.1 3.1 

1.590 -0.750 247.5 5.63 0.419 1.22 29.3 4.8 2.1 

1.590 -0.750 292.5 5.64 0.417 1.28 30.0 4.8 1.5 

1.590 -0.750 337.5 5.62 0.420 1.32 37.9 5.2 2.3 

1.590 -0.583 22.5 5.61 0.422 1.31 45.4 6.8 2.4 

1.590 -0.583 67.5 5.64 0.418 1.28 34.8 6.1 2.9 

1.590 -0.583 112.5 5.64 0.417 1.25 32.8 5.9 3.5 

1.590 -0.583 157.5 5.63 0.418 1.22 38.8 5.5 2.8 

1.590 -0.583 202.5 5.63 0.418 1.22 33.6 5.1 2.8 

1.590 -0.583 247.5 5.63 0.419 1.19 43.3 6.6 3.3 

1.590 -0.583 292.5 5.64 0.418 1.24 39.1 6.2 2.6 

1.590 -0.583 337.5 5.62 0.421 1.34 39.1 6.3 2.3 

1.590 -0.417 22.5 5.62 0.421 1.27 35.5 8.0 3.7 

1.590 -0.417 67.5 5.65 0.417 1.15 16.3 5.9 1.8 

1.590 -0.417 112.5 5.64 0.418 1.09 23.0 6.0 2.7 

1.590 -0.417 157.5 5.63 0.417 1.19 43.1 6.4 3.9 

1.590 -0.417 202.5 5.63 0.419 1.19 37.9 6.0 2.8 

1.590 -0.417 247.5 5.65 0.416 1.08 36.5 7.2 2.3 

1.590 -0.417 292.5 5.65 0.417 1.16 30.9 6.9 1.9 

1.590 -0.417 337.5 5.62 0.421 1.29 33.5 7.5 2.4 

1.590 -0.250 22.5 5.61 0.423 1.26 32.9 9.6 2.4 

1.590 -0.250 67.5 5.66 0.416 1.08 21.6 7.8 2.1 

1.590 -0.250 112.5 5.65 0.416 1.06 22.7 7.1 2.9 

1.590 -0.250 157.5 5.63 0.418 1.18 35.6 7.0 2.1 

1.590 -0.250 202.5 5.62 0.419 1.15 33.2 6.6 2.1 

1.590 -0.250 247.5 5.66 0.415 1.01 31.9 7.8 4.0 

1.590 -0.250 292.5 5.65 0.417 1.11 38.7 8.7 4.1 

1.590 -0.250 337.5 5.61 0.424 1.27 35.7 9.0 2.4 

1.590 -0.083 22.5 5.59 0.427 1.18 21.0 9.2 3.2 

1.590 -0.083 67.5 5.65 0.417 0.95 19.3 9.2 1.2 

1.590 -0.083 112.5 5.68 0.414 0.86 22.3 8.9 2.0 

1.590 -0.083 157.5 5.64 0.416 1.13 32.4 7.7 1.9 

1.590 -0.083 202.5 5.63 0.418 1.12 40.1 7.2 2.9 

1.590 -0.083 247.5 5.67 0.416 0.89 32.8 9.7 4.3 

1.590 -0.083 292.5 5.67 0.415 0.98 27.1 9.3 5.0 

1.590 -0.083 337.5 5.59 0.427 1.21 33.4 9.9 3.1 

1.590 0.083 22.5 5.58 0.431 1.15 29.2 9.9 1.9 

1.590 0.083 112.5 5.64 0.419 0.75 26.2 11.5 6.0 

1.590 0.083 157.5 5.63 0.419 1.12 33.7 8.3 2.3 

1.590 0.083 202.5 5.63 0.417 1.11 35.6 8.1 1.9 

1.590 0.083 247.5 5.68 0.415 0.79 35.4 13.2 4.7 

1.590 0.083 292.5 5.67 0.417 0.88 48.8 14.3 5.7 

1.590 0.083 337.5 5.57 0.431 1.16 55.6 11.9 5.3 

1.590 0.250 22.5 5.58 0.431 1.07 30.3 11.5 4.6 

1.590 0.250 112.5 5.65 0.417 0.66 24.2 13.1 4.1 

1.590 0.250 157.5 5.63 0.417 1.05 25.0 8.9 3.7 

1.590 0.250 202.5 5.63 0.418 1.06 36.7 8.7 3.3 

1.590 0.250 247.5 5.67 0.416 0.70 37.9 13.0 3.8 

1.590 0.250 292.5 5.68 0.415 0.87 22.7 11.0 2.2 



d*cr 

[ — 1 
L sr -I 



y stat ^syst 
r nb 1 r nb i 



1.590 0.250 337.5 5.59 0.429 1.08 41.8 11.4 3.9 

1.590 0.417 22.5 5.59 0.429 1.09 45.2 13.0 3.1 

1.590 0.417 67.5 5.68 0.414 0.77 39.0 14.3 4.0 

1.590 0.417 112.5 5.67 0.415 0.68 39.2 13.8 5.5 

1.590 0.417 157.5 5.65 0.416 1.05 27.4 9.9 2.4 

1.590 0.417 202.5 5.65 0.416 1.07 40.7 9.8 2.5 

1.590 0.417 247.5 5.64 0.419 0.75 25.8 12.8 2.4 

1.590 0.417 292.5 5.67 0.417 0.86 52.8 14.9 5.9 

1.590 0.417 337.5 5.58 0.431 1.06 48.7 12.7 3.7 

1.590 0.583 22.5 5.65 0.419 0.92 51.5 17.1 4.5 

1.590 0.583 67.5 5.71 0.411 0.71 38.5 16.7 4.7 

1.590 0.583 112.5 5.67 0.415 0.67 34.2 14.7 4.9 

1.590 0.583 157.5 5.64 0.417 1.00 38.1 11.8 3.5 

1.590 0.583 202.5 5.64 0.417 1.02 31.6 10.3 2.1 

1.590 0.583 247.5 5.66 0.415 0.76 47.7 14.8 4.0 

1.590 0.583 292.5 5.71 0.410 0.81 52.4 17.0 4.1 

1.590 0.583 337.5 5.65 0.420 0.95 31.3 13.1 5.4 

1.590 0.750 22.5 5.60 0.426 0.92 58.3 22.2 13.0 

1.590 0.750 67.5 5.63 0.420 0.85 29.3 15.9 4.1 

1.590 0.750 112.5 5.65 0.417 0.89 30.4 11.5 4.0 

1.590 0.750 157.5 5.64 0.417 1.05 29.7 11.1 2.8 

1.590 0.750 202.5 5.64 0.417 1.05 35.4 10.1 2.6 

1.590 0.750 247.5 5.65 0.416 0.94 39.1 11.9 3.1 

1.590 0.750 292.5 5.63 0.421 0.92 61.9 16.8 3.5 

1.590 0.750 337.5 5.61 0.426 0.91 67.8 20.9 7.2 

1.590 0.917 67.5 5.66 0.415 1.03 25.0 14.4 4.0 

1.590 0.917 112.5 5.66 0.415 1.03 23.5 12.8 2.4 

1.590 0.917 247.5 5.65 0.417 1.05 24.6 12.1 1.6 

1.590 0.917 337.5 5.64 0.421 1.05 39.8 18.0 3.0 

1.625 -0.917 22.5 5.56 0.414 1.41 33.6 4.7 2.2 

1.625 -0.917 67.5 5.56 0.415 1.36 33.2 4.8 3.3 

1.625 -0.917 112.5 5.55 0.416 1.34 36.3 5.0 2.1 

1.625 -0.917 157.5 5.55 0.415 1.32 24.5 4.4 1.8 

1.625 -0.917 202.5 5.55 0.416 1.30 33.9 4.6 2.2 

1.625 -0.917 247.5 5.55 0.415 1.32 23.4 4.0 1.7 

1.625 -0.917 292.5 5.55 0.416 1.33 32.2 4.4 2.0 

1.625 -0.917 337.5 5.56 0.415 1.40 28.2 4.3 2.8 

1.625 -0.750 22.5 5.54 0.419 1.35 37.2 5.5 2.2 

1.625 -0.750 67.5 5.56 0.414 1.30 27.7 4.8 1.7 

1.625 -0.750 112.5 5.56 0.415 1.24 37.3 5.5 2.9 

1.625 -0.750 157.5 5.55 0.415 1.23 35.4 4.9 2.5 

1.625 -0.750 202.5 5.55 0.415 1.21 24.8 4.1 1.3 

1.625 -0.750 247.5 5.55 0.416 1.21 36.4 4.9 2.6 

1.625 -0.750 292.5 5.57 0.414 1.28 23.5 4.2 1.8 

1.625 -0.750 337.5 5.54 0.419 1.33 30.4 4.8 1.7 

1.625 -0.583 22.5 5.53 0.420 1.33 35.0 6.3 2.2 

1.625 -0.583 67.5 5.57 0.415 1.24 25.6 5.9 3.6 

1.625 -0.583 112.5 5.55 0.417 1.17 23.1 5.4 2.3 

1.625 -0.583 157.5 5.54 0.417 1.20 34.6 4.9 2.3 

1.625 -0.583 202.5 5.55 0.415 1.21 30.9 4.6 1.9 

1.625 -0.583 247.5 5.56 0.415 1.13 18.2 4.7 2.2 

1.625 -0.583 292.5 5.56 0.416 1.21 18.3 4.9 1.3 

1.625 -0.583 337.5 5.54 0.419 1.30 29.5 6.1 2.5 

1.625 -0.417 22.5 5.52 0.423 1.25 22.3 7.7 1.7 

1.625 -0.417 157.5 5.55 0.415 1.14 28.7 5.6 1.7 

1.625 -0.417 202.5 5.55 0.416 1.15 23.3 5.1 1.2 

1.625 -0.417 247.5 5.57 0.415 0.93 23.8 6.7 3.1 

1.625 -0.417 292.5 5.57 0.414 1.06 14.9 6.5 2.1 

1.625 -0.417 337.5 5.52 0.423 1.29 22.7 7.6 1.6 

1.625 -0.250 67.5 5.57 0.417 0.90 16.6 8.9 3.2 

1.625 -0.250 157.5 5.56 0.414 1.15 25.5 6.6 1.4 



29 



W cos8* 
[GeV] 



[deg. 



(QlJ <c> 

rGeV 2 i 



rad . 
corr. 



d*cr 
dS2* 

[ — ] 

L sr -I 



Cstat Osyst 
r nb l r nb i 



1.625 -0.250 202.5 5.56 0.414 1.10 23.1 5.8 1.2 

1.625 -0.250 247.5 5.56 0.416 0.85 38.6 10.8 5.0 

1.625 -0.250 292.5 5.56 0.417 0.96 21.7 9.3 2.5 

1.625 -0.083 157.5 5.57 0.412 1.08 15.5 7.1 1.0 

1.625 -0.083 202.5 5.57 0.413 1.07 14.4 6.5 1.5 

1.625 0.083 157.5 5.56 0.415 1.05 13.5 8.3 1.1 

1.625 0.083 202.5 5.56 0.414 1.07 18.8 7.6 1.4 

1.625 0.250 157.5 5.56 0.414 1.02 15.9 8.9 2.4 

1.625 0.250 202.5 5.56 0.414 1.00 19.2 8.7 1.6 

1.625 0.417 202.5 5.57 0.412 0.98 23.7 9.9 2.0 

1.625 0.583 202.5 5.57 0.413 0.96 20.3 10.7 1.3 

1.625 0.750 112.5 5.62 0.409 0.63 23.3 13.3 1.7 

1.665 -0.917 22.5 5.46 0.413 1.39 24.6 4.6 2.0 

1.665 -0.917 67.5 5.46 0.412 1.38 20.6 4.7 1.5 

1.665 -0.917 112.5 5.45 0.413 1.34 19.3 4.7 2.0 

1.665 -0.917 157.5 5.45 0.413 1.32 26.1 4.6 1.7 

1.665 -0.917 202.5 5.45 0.413 1.33 15.3 4.0 2.0 

1.665 -0.917 247.5 5.45 0.413 1.33 23.9 4.6 1.9 

1.665 -0.917 292.5 5.46 0.412 1.35 31.3 4.9 3.4 

1.665 -0.917 337.5 5.46 0.413 1.37 21.6 4.3 2.0 

1.665 -0.750 22.5 5.45 0.415 1.34 26.5 6.1 1.7 

1.665 -0.750 67.5 5.46 0.412 1.31 14.8 4.7 1.7 

1.665 -0.750 112.5 5.46 0.412 1.23 17.3 4.7 1.5 

1.665 -0.750 157.5 5.46 0.412 1.23 21.8 4.6 1.7 

1.665 -0.750 202.5 5.46 0.412 1.23 22.2 4.3 1.5 

1.665 -0.750 247.5 5.45 0.412 1.20 24.8 4.9 1.9 

1.665 -0.750 292.5 5.46 0.412 1.25 25.6 4.8 1.6 

1.665 -0.750 337.5 5.44 0.416 1.33 22.7 5.5 2.5 

1.665 -0.583 22.5 5.42 0.421 1.31 17.6 7.4 2.8 

1.665 -0.583 67.5 5.48 0.412 1.18 16.2 7.2 3.0 

1.665 -0.583 112.5 5.48 0.410 1.17 12.7 5.7 2.9 

1.665 -0.583 157.5 5.45 0.413 1.21 16.6 5.0 2.0 

1.665 -0.583 202.5 5.45 0.413 1.19 18.4 4.7 1.0 

1.665 -0.583 247.5 5.49 0.408 1.11 14.0 5.9 1.7 

1.665 -0.583 292.5 5.46 0.413 1.19 21.3 7.3 2.0 

1.665 -0.583 337.5 5.42 0.420 1.29 17.2 7.4 1.6 

1.665 -0.417 157.5 5.46 0.411 1.15 15.2 5.6 1.2 

1.665 -0.417 202.5 5.46 0.411 1.16 13.1 5.2 1.1 

1.665 -0.250 157.5 5.47 0.409 1.12 11.4 6.4 2.0 

1.665 -0.250 202.5 5.47 0.410 1.10 10.4 5.9 2.0 

1.705 -0.917 22.5 5.35 0.411 1.39 20.8 5.0 1.4 

1.705 -0.917 67.5 5.36 0.409 1.35 18.4 5.1 1.6 

1.705 -0.917 112.5 5.36 0.409 1.35 21.2 5.5 1.8 

1.705 -0.917 157.5 5.35 0.410 1.33 22.0 4.8 1.5 

1.705 -0.917 202.5 5.35 0.410 1.32 21.2 4.6 1.4 

1.705 -0.917 247.5 5.35 0.410 1.31 26.5 5.3 2.0 

1.705 -0.917 292.5 5.35 0.409 1.33 18.5 4.5 1.2 

1.705 -0.917 337.5 5.35 0.411 1.38 24.6 4.8 1.8 

1.705 -0.750 22.5 5.34 0.412 1.36 22.7 7.4 1.7 

1.705 -0.750 67.5 5.36 0.409 1.27 12.1 5.2 1.6 

1.705 -0.750 112.5 5.36 0.408 1.23 15.9 5.3 2.2 

1.705 -0.750 157.5 5.35 0.409 1.24 20.6 5.0 1.4 

1.705 -0.750 202.5 5.35 0.409 1.24 11.9 4.3 1.5 

1.705 -0.750 247.5 5.35 0.409 1.19 17.2 5.3 2.8 

1.705 -0.750 292.5 5.35 0.410 1.26 23.2 5.4 1.6 

1.705 -0.750 337.5 5.34 0.412 1.35 21.8 6.5 3.1 

1.705 -0.583 157.5 5.35 0.410 1.21 13.8 5.8 1.1 

1.705 -0.583 202.5 5.35 0.410 1.19 11.9 5.4 1.8 

1.745 -0.917 22.5 5.25 0.408 1.41 19.0 5.1 2.1 

1.745 -0.917 67.5 5.25 0.407 1.36 16.6 5.3 2.1 

1.745 -0.917 112.5 5.25 0.407 1.34 8.2 4.7 1.1 



W cos9* 
[GeV] 



[deg 



V (Qbin) <£> 
GcV 2 i 



rad . 
corr. 



d*cr 
dS2* 

[ — ] 

L sr -I 



y stat ^syst 
r nb l r nb "i 



1.745 -0.917 157.5 5.24 0.408 1.32 13.5 4.7 1.3 

1.745 -0.917 202.5 5.24 0.407 1.32 13.3 4.6 1.5 

1.745 -0.917 247.5 5.25 0.407 1.34 18.7 5.0 2.6 

1.745 -0.917 292.5 5.25 0.407 1.37 15.5 4.7 1.2 

1.745 -0.917 337.5 5.24 0.409 1.38 17.4 4.8 1.3 

1.745 -0.750 22.5 5.22 0.412 1.42 20.9 8.6 4.1 

1.745 -0.750 67.5 5.26 0.406 1.29 16.0 5.8 2.2 

1.745 -0.750 112.5 5.26 0.407 1.28 12.7 5.5 1.0 

1.745 -0.750 157.5 5.25 0.406 1.26 14.1 5.0 1.6 

1.745 -0.750 202.5 5.25 0.407 1.24 9.0 4.4 0.7 

1.745 -0.750 247.5 5.25 0.407 1.15 9.2 5.2 0.9 

1.745 -0.750 292.5 5.25 0.408 1.22 13.7 5.7 2.4 

1.745 -0.750 337.5 5.23 0.412 1.33 16.2 8.2 4.7 

1.785 -0.917 22.5 5.13 0.407 1.37 19.2 5.6 4.4 

1.785 -0.917 112.5 5.14 0.405 1.36 11.5 5.4 1.0 

1.785 -0.917 247.5 5.14 0.406 1.31 8.8 4.8 0.7 

1.785 -0.917 337.5 5.13 0.407 1.38 13.5 5.1 1.9 

1.785 -0.750 157.5 5.14 0.404 1.26 14.3 5.4 1.3 

1.830 -0.917 22.5 5.02 0.406 1.37 18.4 5.9 1.7 

1.830 -0.917 67.5 5.03 0.404 1.37 12.5 5.0 1.4 

1.830 -0.917 112.5 5.03 0.404 1.35 17.8 6.2 1.8 

1.830 -0.917 247.5 5.02 0.405 1.29 13.4 5.6 1.3 

1.830 -0.917 292.5 5.02 0.405 1.35 10.4 5.0 1.2 

1.830 -0.917 337.5 5.02 0.405 1.38 9.0 5.1 2.4 

1.830 -0.750 157.5 5.03 0.403 1.27 11.4 5.8 0.8 



TABLE XIII: Higher-Q 2 extracted differential cross-section. 



W cos^; <£„ {QiJ (e) 



[GeV] 



[deg. 



rGcVf 



corr. dn* ° stat 
r nb I r nb 1 



Jsyst 
" nb 1 



1.500 -0.833 36.0 7.04 0.216 1.52 49.4 14.2 3.4 

1.500 -0.833 108.0 7.06 0.211 1.56 32.7 7.6 1.5 

1.500 -0.833 180.0 7.07 0.210 1.57 36.9 7.8 2.0 

1.500 -0.833 252.0 7.06 0.210 1.54 38.9 7.8 2.3 

1.500 -0.833 324.0 7.04 0.216 1.53 36.4 12.3 2.0 

1.500 -0.500 108.0 7.06 0.211 1.53 35.8 8.1 2.0 

1.500 -0.500 180.0 7.07 0.210 1.53 36.5 7.8 1.8 

1.500 -0.500 252.0 7.06 0.211 1.53 42.6 7.9 2.1 

1.500 -0.167 108.0 7.07 0.210 1.52 33.8 9.0 1.9 

1.500 -0.167 180.0 7.08 0.208 1.50 40.8 10.5 2.5 

1.500 -0.167 252.0 7.07 0.209 1.52 50.0 10.0 3.3 

1.500 0.167 180.0 7.07 0.209 1.45 38.0 15.3 2.7 

1.500 0.167 252.0 7.06 0.212 1.48 38.5 13.6 2.2 

1.530 -0.833 36.0 6.97 0.216 1.30 24.2 11.3 2.1 

1.530 -0.833 108.0 7.01 0.209 1.40 38.5 6.6 2.3 

1.530 -0.833 180.0 7.01 0.208 1.37 42.6 6.5 2.3 

1.530 -0.833 252.0 7.01 0.209 1.38 55.2 7.4 2.7 

1.530 -0.833 324.0 6.97 0.216 1.31 43.4 14.2 2.9 

1.530 -0.500 108.0 7.00 0.210 1.38 43.1 7.8 1.9 

1.530 -0.500 180.0 7.01 0.208 1.36 34.8 6.7 1.9 

1.530 -0.500 252.0 7.01 0.210 1.37 31.6 6.3 1.5 

1.530 -0.167 108.0 7.02 0.208 1.38 37.4 8.5 1.8 

1.530 -0.167 180.0 7.05 0.204 1.33 47.1 10.5 2.3 

1.530 -0.167 252.0 7.02 0.208 1.37 28.7 7.3 1.7 

1.560 -0.833 108.0 6.93 0.208 1.27 34.5 6.0 1.6 

1.560 -0.833 180.0 6.94 0.207 1.26 42.4 6.0 2.1 

1.560 -0.833 252.0 6.93 0.208 1.26 37.9 5.8 1.9 
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rad. 
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Osyst 


[GeV] 




[deg.] [*£] 






[ — ] 

L sr J 


[ — ] 

L sr J 


[— 1 

L sr -1 


1.560 


-0.500 


108.0 


6.94 


0.207 


1.25 


35.6 


6.5 


1.7 


1.560 


-0.500 


180.0 


6.95 


0.206 


1.24 


37.5 


6.2 


1.7 


1.560 


-0.500 


252.0 


6.94 


0.208 


1.22 


40.9 


6.6 


1.9 


1.560 


-0.167 


108.0 


6.95 


0.205 


1.22 


21.8 


7.4 


1.2 


1.560 


-0.167 


180.0 


6.96 


0.204 


1.25 


33.0 


10.9 


2.1 


1.560 


-0.167 


252.0 


6.96 


0.205 


1.24 


22.3 


6.8 


1.1 


1.560 


0.167 


252.0 


6.97 


0.203 


1.03 


50.8 


19.5 


9.1 


1.592 


-0.833 


108.0 


6.85 


0.207 


1.18 


19.3 


4.5 


1.3 


1.592 


-0.833 


180.0 


6.86 


0.206 


1.16 


27.5 


4.6 


1.2 


1.592 


-0.833 


252.0 


6.85 


0.207 


1.17 


20.5 


4.3 


0.9 


1.592 


-0.500 


108.0 


6.86 


0.205 


1.19 


19.8 


5.1 


1.0 


1.592 


-0.500 


180.0 


6.87 


0.204 


1.19 


24.1 


5.1 


1.2 


1.592 


-0.500 


252.0 


6.86 


0.205 


1.16 


22.6 


5.7 


1.4 



1.592 


-0.167 


108.0 


6.89 


0.202 


1.05 


29.5 


10.1 


1.5 


1.635 


-0.833 


108.0 


6.74 


0.205 


1.15 


28.0 


5.4 


1.5 


1.635 


-0.833 


180.0 


6.75 


0.203 


1.15 


21.8 


4.3 


1.1 


1.635 


-0.833 


252.0 


6.74 


0.205 


1.17 


16.9 


4.2 


0.9 


1.635 


-0.500 


108.0 


6.76 


0.203 


1.13 


16.5 


5.3 


0.9 


1.635 


-0.500 


180.0 


6.75 


0.203 


1.18 


10.3 


4.6 


0.6 


1.685 


-0.833 


108.0 


6.60 


0.203 


1.26 


11.8 


5.6 


0.7 


1.685 


-0.833 


252.0 


6.60 


0.203 


1.27 


12.0 


4.9 


0.9 


1.685 


-0.500 


180.0 


6.60 


0.202 


1.19 


15.8 


4.9 


0.7 


1.740 


-0.833 


180.0 


6.46 


0.201 


1.30 


18.3 


4.7 


1.0 
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